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PART I. LIDAR CROSS SECTION OF A TILTED CUBE-CORNER RETROREFLECTOR 

ABSTRACT 

Lidar cross section of some typical types of cube-corner retroreflectors 
(CCR's) having a three corner mirror system is investigated for the case that 
the CCR is tilted from the normal illumination axis. Analytic expressions for 
the effective aperture area for the tvo typical window types (circular and 
hexagonal) of CCRs are obtained for the case that the CCR is tilted. The range 
of incidence angle in which only the total reflection occurs at all three 
uncoated corner mirrors has been found to vary considerably with the 
orientation of CCR and the refractive index of the CCR prism. The analytical 
expression for the far-field diffraction pattern of a tilted CCR is obtained by 
taking different polarization transformation of the six sectors of the 
effective reflecting aperture into account. This expression is essential when 
evaluating the lidar cross section of a moving CCR which is tilted in general. 
Formulas for the angles defining the six sectors have also been obtained. 
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I. INTRODUCTION 

Analysis of the lidar cross section of a cube corner retroref lector (CCR) 
has become important under the enlightenment of its use in ranging satellites 
[1] and on the moon [2]. There are various types of CCRs or array of CCRs. It 
may be designed to have a hexagonal, triangular, or circular front view. 
According to the existence of coating on the triple-corner-mirror system, they 
may be classified to coated-mirror CCRs for wider retroreflecting solid angle 
and uncoated-mirror CCRs for better reflectivity and longer life in space. 
According to the array shape of CCRs, there are honey-comb array (retroref lec- 
tors for traffic signs and automobiles), flat array (on the moon), spherical 
array (e.g., the Laser Geodynamic Satellite (LAGEOS)), etc. Some CCRs have an 
intentional dihedral angle error for diffusing the reflected beam. In every 
case, lidar cross section analysis gives us some ideas on how to analyze the 
signal received from existing CCR arrays and on how to design those new CCRs. 

We shall consider the lidar cross section of a single CCR whose front view 
is either circular or hexagonal, whose triple-mirror system is either coated or 
uncoated. To this end, the range of effective angle of incidence and the area 
of effective retroreflecting aperture is found. The far-field diffraction 
effect will also be considered as well as the velocity aberration effect. The 
analytical expression of the far-field diffraction pattern is obtained for a 
tilted CCR in terms of the Kirchhoff integral of the complex amplitudes of two 
linear orthogonal polarizations transformed while being reflected. The 
effective reflecting aperture over which the Kirchhoff integral should be 
carried is divided into six sectors having their own transformation matrices. 

It is a generalization of the work of Chang et al. [2] who obtained the 
expression of the diffraction pattern for a normally illuminated CCR with a 
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circular face. Throughout this work the effect of dihedral angle error is 
assumed to be nonexistent. 

II. EFFECTIVE RETROREFLECTING APERTURE AND REFLECTIVITY 
Suppose a CCR with a circular front window, sometimes called a cylindrical 
CCR, is tilted from the normal illumination. That is, light enters the CCR 
with the entrance angle *. Then, the portion of effective retroreflecting 
aperture looks like Fig. 1 if we look at the CCR along the illumination axis, 
and its area is given by 

A($) = 2a^ [cos - *(^ tan$) - ^ tan<|> 7\ - (-)^tan^<j)J cost, (1) 

3 3 3 . 


9 = sin ( — - — ) 


( 2 ) 


where a is the radius of the circular face and d is the depth of the vertex 
from the front face, n is the refractive index of the CCR prism. For the 
front face to be circular, d > /2a. Usually d = /2a is chosen for maximum 
A ($) and henceforth ve shall assume this configuration for a CCR with a 
circular face. The formula for d = /2a was given in Re f . [ 1 ] . Note that the 
reflectivity and the area of the effective reflecting aperture are important 
factors that determine retroref lectance in the near-field. In fact, if the 
triple-corner-mirror system in a CCR is coated with metal, the reflectivity 
depends on the kind of metal as well as on the angle of incidence and the 
polarization of the incident light. Vhereas, if corner mirrors are uncoated 
the reflectivity depends on whether the condition of the total reflection is 
satisfied or not inside the CCR prism. If total reflection occurs three times 
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(a) 



Figure 1. Two front views and the side view of a cube-corner ret roref lector 
with a circular face, a) the front view when the observation axis 
is coincident with the reference axis which is perpendicular to the 
window face, b) the front view when the observation axis is tilted 
from the reference axis. The effective reflecting aperature is 
shown unhatched with the coordinates (x,y) whose origin is at the 
vertex seen through the window face, c) the side view with the 
image (dashed) of the front window that would be seen inside 
through itself. 
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consecutively inside the CCR prism, the reflectivity is unity during those 
three reflections. Hence the light intensity is only attenuated as light 
transmits into and out of the CCR prism through the front face. However, if 
the condition of the total reflection is not satisfied at one of the three 
reflections, the reflectivity drops to a negligibly small fraction and 
considerable light intensity passes through that surface. Such a negligibly 
small fraction of retroreflection is due to the Fresnel reflection at the 
corner mirror where total reflection does not happen. 

Here we shall consider the condition of the total reflection associated 
with the direction of the illumination axis and the orientation of a CCR. We 
define various directions and their corresponding angles as they are shown in 
Fig. 2. For simplicity, we shall assume that there is no dihedral angle error, 
that is, the three corner mirrors are at right angle to each other. To find 
the condition of the total reflection we have to find the incidence angle of 
light with respect to each of the three corner mirrors where the light is 
reflected consecutively inside the CCR prism. Two angles are given — *, the 
incidence angle of light illumination on the front window, and y, the azimuth 
angle measured from the direction to one of the three azimuths for the three 
edges of a CCR. At each edge two corner mirrors intersect. Given these two 
angles, we can find a and 3> the latitude and longitude, respectively, in the 
spherical coordinate system (Appendix A). 

sina = — (cos $ + /l sin<fr cosy), 0 < a < rt/2, (3) 

/3 

cos(fS _ ", = /3 j ;ost - sina , 

sir/a 


0 < 3 < it/2, 


(4) 
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where <j> is the refracted angle inside the CCR prism given in Eq. (2). We 
may also find cosine of the incidence angles of light on the three corner- 
mirror planes. They are sina for the mirror plane placed horizontally in the 
coordinate system shown in Fig. 2, and cosa cos(3 and cosa sin0 for the two 
mirrors placed vertically and at right angle to each other. When the corner 
mirrors are not coated, all of these three values of cosine should be less than 

J\ - 1/n^, cosine of the critical angle of total reflection, to ensure the 
total reflection. That is, 


cosY^ 
cos ^ 
cosY^ 


■ cosa cos0 < J\ - 1/n^, 


* cosa sing < /l - 1/n^, 
= sina < /l - 1/n^, 


(5) 

( 6 ) 

(7) 


These three inequalities determine the retroreflecting angle range of the 

incidence angle $ at a specific orientation of the CCR. Only one of the three 

incidence angles does not satisfy the condition when the three inequalities are 

not satisfied in all. Let this incidence angle be Y . Then for the TE 

m 

polarization whose electric field is parallel to the mirror plane and the TM 
polarization which is orthogonal to TE, the reflection coefficients for the 
electric field amplitude are 


tan <Y m - Y m ) 

^TM^m^ = tan (Y + Y ) ’ ‘WV 
m m 


sin ( *m - V 

sin (Y m + Y m ) 


( 8 ) 


respectively, where Y = sin ^ (n sin Y ). 

m m 

If we neglect the retroref lection by multiple-order reflection inside the 
CCR prism. The first-order reflectivity of an uncoated-mirror CCR is 
determined by the attentuation when light passes through the front window in 
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and out. This attentuation depends on the polarization of the incident light 
and on the existence of antireflection coating. Here we shall not include the 
effect of antireflection coating in our analysis. For the TE polarization 
whose electric field is parallel to the front window and the TM polarization 
which is orthogonal to TE, the transmission coefficients for the electric field 
amplitude as it enters the CCR are [3] 


in 

t (*) 
TE 


2 cos$ 

cost + ncos<)> * 


in 

t (*) 
TM 


2 cos$ 

ncos$ + cos<(> 


(9) 


respectively. Those coefficients for the light coming out are 


out 

t (*) 
TE 


2n cos4> 

n cos <t> + cost ’ 


out 

t (*) 
TM 


2n cos<j> 
cos<f> + n co si 


( 10 ) 


Here we have not considered the polarization change brought about during 
the three consecutive total reflections at the three corner mirrors. It is 
expected that the retroreflected light will come out of the front window with 
six different polarizations depending on the location where the light comes out 
[ 2 J . One may observe those six sectors that look like pie cuts as he sees 
through a tilted CCR. 

For a given aspect of the CCR, there are three incidence angles, ^ , ^ 2 * 
and ^ inherent to the reflections at the three corner mirrors regardless of 
their order of reflection. However, light comes out with six different 
polarization states according to the six different permutations of the three 
internal reflections. It is just as a series of different birefingence plates 
do not commute unless their corresponding two orthogonal eigenpolarization 
states coincide between the plates [4J. 



If the triple-mirror system is coated with metal, e.g., aluminum or 
silver, the expression for the attenuation in this case will be rather 
complicated than the case of an uncoated CCR, and the result may be inaccurate 
unless very accurate data of the conductivity and the dielectric constant of 
the metal at the wavelength of the light being used is provided. 

Most CCRs that are used as retroreflectors in vehicles or for traffic 
signs have a hexagonal front view and form a flat honeycomb-like array made of 
plastics without rear-side coating. Hence the knowledge on the effective 
reflecting area A($, y) of a tilted CCR with a hexagonal front window is needed 
for the calculation of its retroreflectance. The method for finding the 
effective reflecting area of a tilted CCR is found in Refs. [5] and [6]. Based 
on this method, we can find explicit expressions for its computation as 
follows. We can easily find that there is a hexagonal symmetry for the 
condition of total reflection and thus for the computation of the effective 
retroref lecting area. So we may consider only one of six equivalent sectors 
after transforming the given azimuth angle into an equivalent azimuth y' in 
this sector. By choosing this sector in 0 < y' < n/3, the transformation is 
given by 


y' = min( | y| , \y-2n/3\, | Y -4tl/3|) . (11) 

Then by the rule given in Eqs.(3) and (4) we can find a new set of angles, a' 
and g' , from and y' . Next, according to the value of tan a', we can find the 
effective retroreflecting aperture area normalized to the case of normal 
incidence ($ = 0). (Appendix B) 


10 


A (»,V) 
A(0,r) 


r 2sin2|S' cosa' cota', 
j for 2cos|3' < tana' , 


cos$ 

/3 sin<|> 




sing' (4cosa' - sina' /cos|3' ) , 
for cos0' < tana' < 2cos{3', 


’*> sin|3' cosa' (4 - cos£' cota'), 
for sing' < tana' < cos3'. 


( 12 ) 


Actually tana' < 1//2 is preexcluded from the possible range of the angle a' 
that the first sector can take. In Fig. 3, this effective retroref lecting 
aperture area and the beam reflectivities are plotted in Fig. 3 for two sample 
azimuthal directions of a CCR by varying the entrance angle. They are plotted 
for two types of CCRs; the one having a circular face with the refractive index 
n = 1.455 of quartz and the other having a hexagonal face with n = 1.5 of 
plastics. The beam reflectivity is the absolute square of the product of the 
two kinds of transmission coefficients in Eqs. ( 9 ) — < 10 ) and the reflection 
coefficient in Eq. (8). The last coefficient may not be considered when total 
reflection only occurs in the triple-mirror system. The FORTRAN program for 
computing and plotting such area or the beam retroreflectivity of the foregoing 
two typical types of CCRs is given in Appendix C. 

The near-field retroreflectance of a CCR with an uncoated corner-mirror 
system can be found by multiplying the effective reflecting aperture area and 
beam reflectivity of light passing through the front window and then 
normalizing it with respect to the case of normal incidence. 

As one can easily see, the effective reflecting area does not vary so much 
with the varying azimuth angle. However, it varies considerably with the 
entrance angle. Moreover, when the CCR has uncoated corner mirrors, the 
acceptance angle range in which only the total reflection occurs at the triple- 
corner-mirror system varies considerably according to the orientation of the 
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Figure 3. 
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(a) 


(a) The variations of the normalized effective retroref lecting 
aperture area and the beam reflectivity with the varying entrance 
angle $ for the cube-corner retroref lector with a circular face. 

The refractive index has been set at n = 1.455 of quartz. Two 
sample asimuth angles y = 0° and 90° from one edge direction 
have been chosen, -oo- The common effective retroref lecting 
aperture area for both azimuths. [- 4 — 4 ] The beam reflectivity 

for y = 0° [y = 90°] when the triple mirror system is uncoated. 

The beam reflectivity for both azimuths when the triple-mirror 

system is coated, (b) Those variations for the one having a 
hexagonal face. The refractive index has been set at 1.5. 

[-6— 6] The effective retroref lecting aperture apea for y = 0° 

[y = 90°]. [-4“4] The beam reflectivity for y = 0° [y = 90°] 

when the tripple-mirror system is uncoated. [-4-4^ The beam 
reflectivity for y = 0° [y = 90°] when the tripple-mirror system is 
coated. 
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Figure 3. (continued) 



NEAR FIELD RETROREFLECTANCE 


Figure 4. 
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Variations of the near field retroref lectance for the two types of 
corner-cube retroref lectors considered in Fig. 3 (a) and (b). 

(a) The one with a circular face (n = 1.455) 

■►> [-A — ArJ for y = 0° [y = 90° ] when uncoated 
*0*0- for both azimuths when coated 

(b) The one with a hexagonal face (n = 1.5) 

> ► [-A-±] for y = Cr [y = 90° ] when uncoated 
for y = 0° [y = 90° ] when coated. 
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Figure 4. (continued) 
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CCR and the refractive index of the CCR prism. In Fig. 5, the maximum 
allowable entrance angle for the three cases of retroref lection by the two 
front-view types of a CCR is plotted as we vary the azimuth angle y. The 
FORTRAN program that computes the maximum entrance angle is the program 
CRITANG, and CRITGRM plots it. They are shown in Appendix C. 

III. LIDAR CROSS SECTION IN THE FAR-FIELD WITH THE VELOCITY ABBERATION EFFECT 

If the CCR is far away from the receiver, we have to consider the 
diffraction of the reflected beam. For a rough estimate, when a circular 
aperture of diameter 2.54 cm is illuminated perpendicularly by light of 
wavelength 6943 A, the first zero in the far-field diffraction pattern occurs 
at 6 arcseconds off-axis from the reference axis. The broader the diffracted 
beam, the less peak power we can detect at a receiving point* If the CCR is 
tilted with respect to the illumination axis, the diffraction pattern becomes 
broader with the less receiving power. A small dihedral angle error may make 
the diffraction pattern irregular, which also decreases the receiving power. 

If the CCR target has a velocity to the direction right to the illumination 
axis, finite round-trip time induces the velocity aberration effect that 
displaces the maximum of the diffraction pattern from the illumination axis. 

In this section, we shall investigate the diffraction pattern of a tilted 
CCR having a circular window, since most CCR arrays being used in space are of 
this type. The diffraction analysis of a tilted CCR has its worth when we 
consider, at the same time, the decrease in receiving power due to the velocity 
aberration effect of the moving target, i.e., the ranging satellite. The 
far-field diffraction pattern for the case of normal incidence to a CCR with a 
circular front window was investigated by Chang et al. [2]. They were 
concerned with the flat rectangular array of CCRs such as the array on the 
moon, which always faces to the earth. 


CRIT. ENTR. ANGLE 



AZIMUTH ANGLE 


Variation of the maximum entrance angle under which the incident 
light is effectively or totally reflected at the triple-mirror 
system as a function of the azimuth angle y- Line (a) and (b) 
represent these variations for the cases of the coated and uncoated 
triple-square-mirror systems, respectively, while line (c) 
represents that for the case of the uncoated triple-mirror system 
with a circular face. 
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If a CCR is tilted, the effective reflecting aperture has a symmetric 

streamlined shape as is shown unhatched in Fig. 1. Suppose the aperture is 

illuminated with a plane wave of the wave number k and the polarization 

represented by its electric field, [E^ x , e^^ z ^ ct ^ where c is the 

velocity of light, and the superscript "T" refers to the transpose. Ve may 

choose the x-component to be the TM component to the front window of the 

CCR, while E^ to be TE component. Then, the far-field diffraction pattern of 

the retroref lec ted light may be found by the Fraunhoffer diffraction formula 

under the Kirchhoff's physical approximation. Considering the six different 

polarization outputs, one may express the electric field at the photodetector 

_2 W Z+ct ) 

on the ground as E (X,Y)e with the column vector, 

o 


E g (X,Y) 


k 

2i 


U 

aperture 


T „ M(9) T. 
out v ' in 


-E. 

IX 

^E. 

iy 


jexp[-i £ (Xx + Yy)] dxdy , (13) 


x = r cos 9 , 


r sin 9 , 


(14) 


aside from the constant phase factor. Here i = T , and L is the distance from 
the CCR to the receiver. (x,y) or (r,9) is the two-dimensional coordinate on 
the CCR aperture, while (X,Y) is that on the receiving plane. Note that the 
x-axis joins the center of the aperture and the center of the front window of 
the CCR on the projection plane. The origin on the (x,y) plane is chosen at 
the center of the effective reflecting aperture. 


T. 


m 


’•TM 

0 



out 


out 

iM 

0 


out 

iE 


(15) 
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Note that there are six different orders of reflection at the three corner 
mirrors, which gives the following six different matrix representations for the 
role of the triple corner-mirror system on polarization according to six 
different ranges of azimuth 0. 


M(9) 


• m 3 m 2 m x , 
m 3 m 2 , 

^ Mj M 3 M 2 , 

M ! M 2 M 3 ’ 

m 2 m x m 3 , 

M 2 M 3 M 1 ’ 


© x + n < 0 < M 3 
© 3 < 0 < © 2 + n 
© 2 + n < 0 < © 3 
© 3 < 0 < © 3 + n 
© 3 + n < 0 < © 2 

© 2 < 0 < 0^ + IT 


(16) 


M. 

J 


f p TM ( ' , 'j ) 
R <- 0 j> ' 0 


«r E <V ] R ^ 


j = 1, 2, 3, 


[ cos0 sin0 , 

. 

-sin© cos0 J 




7 2 

cos\I/. - in /n sin2i{;. - 1 
J 3 


costJ / 


*TE ( V = 


ncos\J/ 


72 2 

. + in /n sin tf/. - 1 
3 3 

2 

- i /n sin - 1 


ncos^ 


. + i v^i^sin^i)/. - 1 


(17) 

(18) 

(19) 


(20) 


where the three angles of incidence to the three corner mirrors, ^ and 
are given as their cosine values in Eqs.(5)-(7), respectively, while the 
angles a and 3 have been found in Eqs.(3)-(4). Various transmission 
coefficients in Eq.(14) are given in Eqs. (9)-(10) . The six sectors of the 
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effective retroreflecting aperture are divided at 9 = 0. , 0. + it, j = 1, 2, 3, 
whose tangents are given by (see Appendix D) 

tan©. = tan©. (cos<|>/cos#) , j = 1 , 2 , 3 . ( 21 ) 

J J 

where ©ys are the azimuth angles dividing the six sectors viewed inside the 
CCR prism as shown in Fig. 6 . They may be found from 


cos© - 1 ~ cosa cosptsina + cos« (cos© + sing)) 

1 2 2 
/3 sin<|> /I - cos a cos p 


( 22 ) 


. 2 2 

a sin a - sina cosa sing + cos a cosg (cosg - sing) 

COS Drt 5 . - ■ _ _ _ __ 

i 1 T 

JZ sin$ /I - cos a sin g 


(23) 


COS 0 _ cosa - sina (cosg + sing) (24) 

^ /3sin$ 

The ambiguity involved in the inverse trigonometry in Eq.(21) can be resolved 
by noting that the angles ©j and ©^ are in the same quadrant in the x-y plane. 

The integration in Eq.(13) may be evaluated by the two-dimensional 
discrete Fourier transform. It may also be benefited by making use of the fact 
that the geometry has a radial symmetry. 

In a special case of a CCR with the three corner mirrors that can preserve 
the polarization as in the case of highly reflecting metal-coated corner 
mirrors, M t (9) can be considered as an identity matrix. In this case, the 
above integration (13) can be evaluated by the extended use of the 
two-dimensional Fourier transform of an isosceles and that of a trapezoid which 
were given in Ref. [7]. That is, we first segment the aperture into 2N stripes 
with an identical width as shown in Fig. 7. The two segments at the top and at 
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Figure 6. The front view of a cube-corner retroref lector seen inside along 

the internal axis (p) refracted from the external observation axis. 

The six sections are defined by the azimuths 9=0,9 + ji, m = 1, 

o t J nr m ’ ’ 



Figure 7. The streamline-shaped effective reflecting aperture is segmented 
into two isosceleses and a number of trapezoids symmetrically 
disposed with respect to the origin. Two outer lines are 
represented by x = + F(y). 


22 


the bottom may be approximated into two isosceles located symmetrically with 
respect to the origin. All the other "2N-2" stripes may be approximated into 
trapezoids. They may be paired symmetrically with respect to the origin so 
that their superposed two-dimensional Fourier transform may be a real function 
on (X,Y). The two outer lines are represented by x = + F(y), where 

F(y) - cos$ (Si ~ y^ - /2 tan<f>). (25) 

Then the integration in (13) aside from premultiplying factors may be expressed 
as a superposition of I^ f n = 1, . . . , N, each of which is the Fourier 
integral of a pair of trapezoids of an identical height h (an isosceles may be 
viewed as a trapezoid with a constricted upper latus). 


P F(y) ; 

J dy / dx exp[-i(£x + Tny)] = 2 



(26) 


where p 


|F _1 <0) | 


-K 


2 tan^<f>, and h = p/N in consequence 


I 

n 


u , n ) = 


sink oo = 


'-g- {sin(x n £; + y n ti) sink[h(h + a n O] 

+ sin(x £ - y h) sink(h(h - a £)]}, 

j n n it 

^4h{x n cos(ynti) sink(hh) 


+ a n y n sink(2y n h) [(cosh(r|/2) - sink(hri)]}, 


sin(oo/2) 

( r 


w t 0 
oo = 0 


1*0 

(27) 


1 = 0 . 


(28) 
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x n = [ F(nh) + F((n-l)h) ]/2, y n = (n - l/2)h (29) 

a n = [F(nh) - F((n-l)h)]/2 . (30) 

The FORTRAN program that computes Eqs. (26)-(30) was given in Appendix E, 
in which these equations were used for finding the lidar cross section of 
each CCR on the LAGEOS. 

If a CCR is moving with the velocity component v' which is perpendicular 
to the illumination axis, the Bradley or velocity aberration effect causes the 
reflected beam pattern to be angularly displaced by an amount [1] 

* v = 2v'/c (31) 

Here the aspect of CCR with respect to the illumination axis is relatively 
unchanged compared to the velocity of the CCR in moving. The lidar cross 
section is defined by [8] 

_ 4 n I(at the receiver) (32) 

11 H( incident irradiance) 

where H and I are the irradiances measured at the CCR and the receiver, 
respectively. 

So far we have assumed that the CCR has no dihedral angle error and no 
other manufacturing errors. In fact, many CCRs have intentional dihedral angle 
errors so that the retroreflected beam may have a little degree of beam spread 
for practical purpose. The CCRs used in recent ranging satellites also have 
some degree of nominal dihedral angle error in the order of one arcsecond to 
have spreaded diffraction pattern at the receiver which is at the transmitter 
site [9]. 
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The spread of the retroreflected beam induced by the dihedral angle error 
S was found to be 3.26 n5, where n is the refractive index of the CCR prism 
[5]. More elaborate work to find the explicit expression for the six diverging 
reflected beams is found in Refs. [10] and [11]. The result of Chandler [11] 
can be incorporated into our analysis of the far field diffraction pattern. 

Effects of some other manufacturing error which is random in nature was 
treated simply by supposing an equation of the form which provides the best fit 
to experimental datas [8]. 

IV. CONCLUSIONS 

In this report, we have analyzed the lidar cross section of one CCR with 
either a circular front window or a hexagonal front window. Analytic 
expressions for the effective reflecting aperture of a CCR with a hexagonal 
window have been found. The polarization change as light is reflected and 
comes out is discussed with a CCR having either uncoated or coated corner 
mirrors. That change may be treated with the Jones calculus. The result of 
this report may directly be in use for CCR arrays either for traffic signs or 
for space use. The expression for the far-field diffraction pattern of a CCR 
has been obtained under the Kirchhoff's physical approximation on aperture 
diffraction, and analysis on the actual pattern to be obtained by the 
polarization analysis is reserved for future research. Besides analysis on the 
effect of the dihedral-angle error has not been included in this report. For 
an array of CCR, the intensity pattern in time for each CCR should be summed 
because of incoherency between reflections by individual CCRs. The far-field 
Fraunhoffer diffraction pattern will be useful in analyzing the lidar cross 
section of an array of CCRs in moving. 


This work has been supported by NASA Goddard Space Flight Center. 
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APPENDIX A: DERIVATION OF EQUATIONS (3)-(4) 


In Figure 2, the reference axis is represented by a unit vector, 

* T 

c = [1//3, 1//3, 1//3] , where the superscript "T" refers to the transpose. 

We shall denote by a unit vector p the direction of light being deflected from 
the illumination axis while entering the CCR window, 


T 

p = [cosa cos|3, cosa sing, sina] 


Then from the cosine rule, c • p = cos$ 


cos <|> = [sina + cosa (cos(3 + sinfl)]//3 


* " rn 

and from |p x c| cosy = p x c • [-1//2, 1//2, 0] , 


(A. 1) 


(A. 2) 


tany = — CQS0C ( cos(3 ~ sin P) 

T 2 sina - cosa (cosS + sing) 


(A. 3) 


From these two equations (A. 2) and (A. 3) we may find Eqs.(3)-(4). 
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APPENDIX B: DERIVATION OF EQUATION(12) 


Assuming no dihedral angle error, one can see, including their images, 
three square mirrors each of which looks like bisecting the other two mirrors 
in the middle. The vertex of the CCR prism is at the center of each of the 
three square mirrors. Define a new two-dimensional Cartesian coordinate system 
for projection whose xy-plane is at right angle to the vector p introduced in 
Appendix A. The vertex of the CCR prism is also chosen as the origin of the 
new coordinate system. One corner along the z-axis would lie on the x-axis on 
the plane of projection. One may assign four coordinates to four points of 
each of the three projected squares. It can be done systematically by the 
following matrix operation: 


,] • e;j] v-»'> r z<«'> 


(B.l) 


T T 

where a is one of eight position vectors such as [Oil] , [01 -1] , etc. 

R^ and R^ describe the passive rotations with respect to the y and z axes, 
respectively. Once all twelve projected coordinates are obtained, one can 
readily find the triply overlapped hexagon and find its area according to one 
of the four conditional phases of the hexagonal shape characterized by the 
corresponding four ranges of tana as were given in association with 
Eq.(12). Finally, considering the index difference, one should multiply 
it with cos$/cos<]>. 
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APPENDIX C 


M.170F 

1 : 

2: 

3: 


6 : 
7: 
3: 
9 : 
10: 
11 : 
12 : 
13: 
14: 
15: 
16: 
17: 
18: 
19: 
20 : 
21 : 
22 : 
23: 
24: 
25: 


2S: 
29: 
30: 
31: 
32: 
33: 
34: 
35: 
36 : 
37: 
38: 
39: 
40: 
41: 
42: 
43: 
44: 
45 : 
46: 
47: 
43: 
49: 
50: 
51: 
52: 
53: 
54: 
5s : 
56: 
57: 
53: 
59: 
60: 
61: 
62: 
63: 
64: 
65: 
66 : 
67: 
68 : 
69: 
70: 
71: 
72: 
73: 


NAME CCRPLT 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c- 


301 


319 


30 

31 
2 

302 



UPDATED ON FEB. 23, 1987 G. HUGH SONG 
makes a graph o' - 

as a function o_ 
with fixed values 

RI: Refractive index of the CCR. 

CCR(DPHI): The resulting relative retro-reflectance normalized 
to the case of normal incidence 
REF : retro-reflectance 
DGAM: See statement No. 302 
DPOL: See statement No. 308 

This program assumes unpolarized input. Should be modified! 


DIMENSION RANGE! 4) ,Y( 180-1) 

DIMENSION APHK180) -AREN ( 180) , CCR (180) 

DATA RANGE/180. ,180. .0. ,1.1/ 

RAD!DE)*DE*ATAN(1. )/45. 

DEG ( RA ) =RA*45 . / ATAN ( 1 . ) 

PI =4 . *ATAN ( 1 . ) 

INPUT 

WRITE! 3,301 ) 

FORMAT! / jJ," ENTER THE FOLLOWING DATA >>>>",/,/,/, 
+ " RI=REFRACTIVE INDEX OF THE CCR.,”,/, 

+ " (THAT OF FUSED SILICA IS 1.455)") 

READ! 3 ,*)RI 


i / / / , 

it an array of hexagonal CCRs (if so, type 1!),".//, 
sincjle CCR with a circular face (if so, type 0!) ? ,//) 


WRITE! 3,319) 

FORMAT!///, 

+ " Is it an arra 
+ " or a 
READ (3 

WRITE! 3,)” NOW, IS THE CCR COATED? YES-1, NO-0” 

READ! 3 , ) LC0AT 
IF ( LCOAt ) 3 0 , 2 , 3 1 

WRITE! 3, ) "WE ASSUME THE METALLIC COAT IS A PERFECT REFLECTOR." 
WRITE! 3, 302) 

FORMAT!/." DGAM ( deg )= ANGLE OF THE PLANE OF INCIDENCE",/, 

+ " MEASURED FROM ONE EDGE DIRECTION (NORTH pole).",/, 

+ " 0. = < DGAM < 180.") 

READ! 3, ) DGAM 
GAM® RAD ( DGAM ) 


3 06 
202 


203 
308 

C 

c 

r* 

c 

c 

C204 

204 
C 


WRITE! 3,306) 

FORMAT! 1,1." IS THE INCIDENT LIGHT POLARIZED? 

READ! 3.) L^OL 

IF ( LP0L >202,204,203 


YES- 




N0-0" ) 


WRITE (3, 308) 

FORMAT!" Sorry. Not developed for polarized light."/, 

+ " The program assumes unpoxarized light in the. following.") 
+" WE SHALL ASSUME THAT THE LIGHT IS LINEARLY POLARIZED . , / , 

+ " SPECIFY THE DIRECTION OF THE POLARIZATION IN deg,",/, 

+ " so that DFOL=90. and DPQL®0. mean the TE and TM incidence. 
+ ,/#" respectively.") 

$eAD( 3 , ) DPOL 

IF ( LPOf. . EQ . 0 ) DP0L=45 . 

DPQL=45 . 


WRITE! 3 ,307 ) 

307 FORMAT!////," INCI. " ,T9 , "Ref ' tv" ,T20 , "REF' ANCE" ,T40 , "EFFECTIVE 

+AREA" , / , " ANGLE" ,T40 , "NORMALIZED" , / / / ) 

DO 81 1=1,180 

APHI ! I ) = 9 0 . -FLOAT ( I ) 

DPH1 = ABS ( APHI ( I ) ) 

IFtI.EQ.90) GAM=GAM+PI 

CALL SUBCCR (RI, GAM, RAD ( DPHI ) ,LFACE, LCOAT, 1, RAD (DPOL) , ALP, 

+EET, DPHIM2 , AREN ( I ) , TRY, CCR! I ) ,T0P) 


28 


/ ** 

75 

76 

77 
73 
79 

ao 

81 

32 

33 
84 
35 
86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 


Us 

117 

118 

119 

120 
121 

123 

124 

125 

126 
127 
123 
123 

130 

131 
122 

133 

134 

135 

136 

137 
133 

139 

140 

141 

142 

143 

144 

145 

146 

147 


81 

7 


83 

84 


100 

101 

102 

103 

C 

104 

105 

106 
107 

5 

108 

109 

110 

C 

111 

997 

112 

113 

114 

995 

115 



C 


IF ( I . EO . 2 ) DPHIM1 =DPHIM2 
DALP=DEG( ALP) 

DBET=DEG( BET) 

WRITE( 36 , 7 ) APHI ( I ) , TRY, CCR ( I) ,AREN( I ) .DALP ,D3ET,CCR( I ) *AREN( I ) 
WRITE ( 3,7 ) APHI ( I ) ,TRY,CCR(I) ,AREN( I) , DALP ,DBET , CCR ( I ) ’“'AREN ( I ) 
CONTINUE 

FORMAT ( 1X,F6 . 2 ,F10 . 6 ,T20 ,F10 . 6 ,T35 ,F10 . 6 ,2 ( 2X,F8 . 3 ) ,2X,F10.6) 


; output and plot 

WRITE ( 3 , 3 0 9 ) DGAM , DPHIM1 , DPHIM2 

309 EFFECTIVE REFLECTANCE OF CCR %%%%%**%*“,//. 

+ //," WHEN THE INCIDENCE PLANE ( DGAM= " , FI 0 . 6 , " deg ) IS FIXED. %/ 
+//'" ACCEPTIBLE ANGLE OF INCIDENCE IS FR0M",F9.6,” TO -",F9.6) 

* 

WRITE( 3,347) 

— fi T.TL 


347 FORMAT ( 


READf 3 , ) 


WHICH TO PLOT? 
NY 


T? Type 1 for the beam relectivity CCR",/ 
OR Type 0 for the norm. eff. area AREN. ) 


IF(NY.EQ.O) 

DO 83- 1 = 1, 130 
Y(I,1)=AREN(I) 
ELSE 

DO 84 1*1,180 
Y(I.1)*CCR(I) 
END tF 


CALL USPLO ( APHI .Y. 180 , 180 , 1,1 , 11HREFLECTANCE , 

+ 11 , 5HANGLE, 5 , 3&REF , 3 , RANGE , 1H1 , 0 , lER 


WRITE (3,) "DO YOU WANT TO RUN ONCE MORE? YES-1, NO-O 
P.EAD( 3, )NY 

IF(NY.EO.O) GO TO 997 

WRITE( 37 ) ” With the same CCR? Yes-1, No-O" 

READ( 3 , ) NY 
IF ( NY ) 5,1,2 


u 

C 

c 


1 

310 


WRITE( 3, ) "WANT A PLOT ON THE ELEC. STATIC PLOTTER? Yes-1 No-O" 

READ (3.) NY 

IF(NY.EQ.l) CALL CCRGRM 

STOP 

END 


SUBROUTINE CCRGRM 

same as the PROGRAM CCRGRM. 

DIMENSION Y ( 180 ) ,RETNEAR( 4 , 130 ) ,APHI(13G) ,CCR(13Q) ,AREN(130) 
REWIND 36 

CALL PINIT( 1.5, 5., 720) 

WRITE( 3.21) 

FORMATC How many sets of curves do you want to plot? (N<5)",/ 
+ ," Each set = norm. eff. retroref. area 8 beam reflty."//, 

+" And you'll plot the near field ref ' ance . " , / ) 

READ (3,) N 
WRITE ( 3,22) 

FORMAT ( s 1st curve : thin solid 2nd : bold solid",/, 

+ " 3rd : thin dahed 4th : bold dahed.") 

DO 3 NPP=1,N 
DO 1 1=1,180 

READ ( 36 . 3 1 0 ) APHI ( I ) , CCR ( I ) , AREN ( I ) , DALP , DBET . RETNEAR ( NFP , I ) 
FORMATC lX,F6 . 2 ,T20 ,1-TO . 6 ,T^5 ,F10 . 6 ,2 ( 2X,F8 . 3 ) ,2X,F10 .6 ) 

IF ( NPP . EO . 1 ) 


CALL MGRPHL( APHI,AP.EN.180,-90. ,90. , ^ . . * . , 
CALL DATPL2( APHI, CCR, 180, 0,1, 0,0. ,0,1,0) 
ELSE IF( (NPP.EO.2) .OR. (NP^.GT.4) ) 

CALL DATPL2 ( APHI , AREN. 180 , 0 , 1 , 0 , 0 . , 0 , 4 , 0 ) 
CALL DATPL2( APHI, CCR, 1 80,6,1 ,6,6. ,6,4,0) 
ELSE IF ( NPF.EQ . 3 ) 

CALL DATPL2 ( APHI , AREN, 180 ,0,1, 0,0. ,0,1,3) 


6. ,4. ) 


29 


CALL DATPL2( APHI, AREN, 18 0,0, 1,0,0. , 0,4,3) 

^CALL^DATPLZ ( APHI , CCR ,180,6,1,6,6. , 0,4^3$ 

SontJnue 

CALL TITLE ( 1HB,21HINCIDENCE ANGLE (deg) , 21 , . 14 , 1 ) 

CALL WHERE 

CALL PINITC 1.5, 5., 720) 

DO 4 ,NPP=1 ,N 
DO 5 1=1,180 
Y(I) =RETNEAR ( NPP , I ) 

IF(NPP.EO.l) CALf, MGRFHLf APHI , Y, 180 , -90 . . 90 . , 0 . , 1 . ,6 . ,4 . ) 
IF ( NPP.EO . 2 ) CALL DATPL2 ( APHI , Y, 180 , 0 ,1 , 6 , 0 . , 0 ,4 , 0 ) 

IF ( NPP . EO . 3 ) CALL DATPL2 ( APHI , Y, 180 ,0 ,1 , 0 ,0 . , 0 , 1 ,3 ) 

IF ( NPP . EQ . 4 ) CALL DATPL2 ( APHI , Y, 180 ,0 , 1 , 0 ,0 . , 0 , 4 , 3 ) 
CONTINUE 

CALL TITLE ( 1HL , 27HNEAR FIELD RETROREFLECTANCE , 27 , . 14 , 1 ) 
CALL TITLE ( 1HB,21HINCIDENCE ANGLE (deg) ,21, . 14,1) 

CALL WHERE 

RETURN 

END 


SUBROUTINE SUBCCR ( RI , GAM , PHI , LFACE , LCOAT , LFRESN ,POL , 

+ ALP, BET, DPHIM,ARN,REF,REFFR ,TQP ) 

G. HUGH SONG, FEB. 21, 1987 

COMPUTES THE REFLECTANCE OF A CUBE CORNER RETRCREFLECTOR 


179: C FROM SPECIFIED DATA OF THE RI, GAM, PHI ANGLES. 

130: C 


30 


>1,92? : 
2 : 
3: 
4: 
5: 
6 : 
7: 
8 : 
9: 
10 : 
11 : 
12: 
13: 
14: 
15: 
16: 
17: 
13: 
19: 
20 : 
21 : 
22 : 
23: 
24: 
25: 


C- 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c- 


FROGRAM CRITANG 


UPDATED ON NOV. 12 f 1986 G. HUGH SONG 
outputs the critical entrance angle,, , 
that depends on the azimuth angle DGAM 
with fixed values RI and DPOL. 

RI: Refractive index of the CCR. , . , 

CCR(DFHI): The resulting relative retro-reflectance normalized 
to the case of normal incidence 
REF : retro-reflectance , 

DGAM: See statement No. 302 
DPOL: See statement No. 303 

SUBROUTINE : SUBCCR 


27: C 

23 : 

29: 

30: 

31: 

32 : 

33: 

34: 

35: 

36: 

37: 

38: 

39: 

40: 

41: 

42: 

43: 

44: 

45: 

46: 

47: 

43: 

49: 

50: 

51: 

52: 

53: 

54: 

55: 

56: 

57: 

58: 

59: 

60: 

61: 

62: 

63: 

64: 

65 : 

66 : 

67: 

68 : 

69 : 

70 : 

71: 


RAD { DE ) =DE*ATAN ( 1 . ) / 45 . 

DEG ( RA ) =RA*45 . /ATAN ( 1 . ) 

PI =4. *ATAN( 1 . ) 

INPUT 

WRi’T'E ( 3 301) 

FORMAT (?./," ENTER THE FOLLOWING DATA >>>>",/,/,/, 
+ " ri=r£f£active INDEX OF THE CCR.,",/. 

+ " (THAT OF FUSED SILICA IS 1.455) r ’) 

P.EAD( 3 , * ) RI 


301 


319 


30 

31 


WRITE ( 3,319) 

FORMAT {III, ^ ^ , __ , . _ . 

" Is it an array of hexagonal CCRs (if so, type l\), ,//, 

+" or a sinale CCR with a circular face (if so, type 01) 7 ,11) 
READ( 3 , * ) LFACE 

WRITE( 3 , ) " NOW, IS THE CCR COATED? YES-1, NC-0" 

READ( 3, ) LCQAT ^ 

5®ITE?3?>" WE 3 ASSUME THE METALLIC COAT IS A PERFECT REFLECTOR." 
GO TO 201 


32 

311 


201 

306 

202 


203 

308 


WRITE; 3,311) „ 

FORMAT ( " When the condition of total reflection is not . /, 

+" satisfied inside the CCR, we ignore the reflectance ,/) 

306) 

FORMAT (/,/," IS THE INCIDENT LIGHT POLARIZED? YES-1, NO-0") 

READ (3.) LPOL 

IF ( LPOf. )202,204,203 

WRITE ( 3,308) 

^CRMAT ( / 

+ ” WE SHAjLl ASSUME THAT THE LIGHT IS LINEARLY POLARIZED.",/, 

+" SPECIFY THE DIRECTION OF THE POLARIZATION IN deg,",/, 

+ •• so that DPOL=90 . and DPOL=0. mean the TE and TM incidence," 
+ , /." respectively.") 

&eAD(3,)DP0L 


204 WRITEf 3,307) _____ 

307 FORMAT (////," AZIMUTH " ,T20 , "CRITICAL ANGLE" ,T40 , "EFFECTIVE 
+AREA" ,111) 

DO 81 IG=0 ,120 1 FORTRAN 77 

DGAM=FLOAT ( IG ) 

GAM=RAD ( DGAM ) 

DELT=45 . 

DFHI=45 

41 CALL SUBCCR (RI, GAM, RAD ( DFHI ), LFACE, LC0AT,0, POL, ALP, 


73 : 


+BET, DPHIM2, AREN, TOTREF, REF .TOP) 
CCR*AR0*ARES*REF* ! Fraunhoffer diff 
IF ( REF . GT . 0 . ) DELT=AES(DELT/2. ) 
IFtREF.LE. 0 . ) DELT=-ABS(DELT/2. ) 

DFHI =DPHI+DELT 

IF( ABS ( DELT) .GT. 1 .E-3 )G0 TO 41 


raction. 





31 


DALP=DEG ( ALP ) 

DBET=DEG ( BET ) 

WRITE (36,310) DGAM , AREN , DPHI , DALP , DBET 
MR ITE ( 3 - 3 1 0 ) DGAM -AREN , DPHI . DALP , DBET 
F0RMAT(1X,F6.2,T20,F10!6,T46,F10.6,2(2X,FS.3) ) 


WRITE(3.) "DO YOU WANT TO RUN ONCE MORE? YES-1, NO-O." 

READ( 3 , )NY 

IF(NY)999, 999,1 

STOP 

END 


SUBROUTINE SUBCCR(RI , GAM , PHI , LFACE ,LC0AT , LFRESN , POL , 
+ ALP, BET, DPHIM,ARN,REF,REFFR,TOP) 


G. HUGH SONG, FEB. 21, 1987 


1 

2 

3 

4 

3 

6 

7 

3 

9 

10 

11 

12 

i 

14 

15 

16 
17 
13 
19 
2C 
21 
22 

23 

24 

25 

2b 

27 



PROGRAM CRITGRM 

DIMENSION X( 0:120) ,Y( 0:120) ,DGAM( 0:120) ,CRIT( 0:120) 
CALL PINIT( 1.5, 5., 600) 

N=0 

DO 1 1*0,120 

READ ( 3 6, 3 1 0 ) DGAM ( I ) , SOMETH , CHIT ( I ) 

CONTINUE 

FORMAT ( 1X,F6 . 2 ,T20 ,F10 . 6 ,T40 ,F1Q . 6 , 2 ( 2X,F8 . 3 ) ) 

DO 2 1=0,120 
X( I) =DGAM( I) 

Y( I ) =CRIT( I ) 

CONTINUE 

IF(N.EQ.O) 

CALL MGRPHL(X,Y, 121, 0.,120., 0.,90., 6., 4.) 

ELSE 

CALL DATPL2 'X,Y,121 , 0,1, 0,0., 0, 3, 0) 

END IF 
N=N+1 

IF ( N . LT . 3 ) GO TO 3 

CALL TITLE ( 1HL, 17HCRIT. ENTR. ANGLE, 17, .14, 1) 
CALL TITLE ( 1HB , 1 3HAZ IMUTH ANGLE, 13, .14, 1) 

CALL WHERE 

STOP 

END 


APPENDIX D: DERIVATION OF EQUATIONS (21)-(24) 


For the moment, we suppose a CCR with unit refractive index. Then we can 
see a tilted front view of a CCR along the observation axis which is the 
illumination axis at the same time. Then we then have the six sectors and six 
angles 9^, 9^ + ft, 92 9^ + n, 9^, and 9^ + it that define the six sectors as 
shown in Fig. 5. Using the cosine rule, we can find the cosine of 9^, 92» and 
9^ as 


c x p • e. x p 

cos9. = ^ 

J ... 

|c x p| le^ x p{ 


j = 1, 2, 3 , 


which will give Eqs . (21)-(23) . Here e^, and e^ are the unit vectors 
directing along the three corners inside a CCR. 

Next, for a CCR with the refractive index n, those 9j's will appear at 
in the ellipce whose short axis is the more shortened by the ratio cos$/cos<j> 


while the long axis remains the same. Therefore 0. and 9^ are related by 
Eq. (21) . 
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> 1 . 7 IP 


APPENDIX E 

NAME DIFPAT 

c 

C G.HUGH SONG , updated on APRIL 15, 1986 

C dives the diffraction pattern for a symmetric pupil function 

C F(y) , which is givn in a function subprogram 

DOUBLE PRECISIONS TOP 

DIMENSION RANGE ( 4 ) ,SPFR( 100) , A (100,1) 

COMMON / C0MF / PHI , TAP 
CHARACTER* 1 XORY 

DATA RI/ 1 . 455/, RANGE/0. ,0. ,0. ,0./ 

RAD ( DEG ) =DEG*ATAN ( 1 . ) / 45 . 

C 

9 WRITE! 3, ) "INPUT DPHI IN deg 1 DPHI=" 

READ ( 3 , ) DPHI 

PHI =RAD( DPHI ) 

SIP=SIN( PHI ) /RI 
TA?=SIP/SQRT( 1 . -SIP*SIP) 

TOP- SORT ( X . -2 4 *TAP*TAP) 

WRITE? 3 , ) " * HCW MANY STRIPES DO YOU TAKE FROM ZERO TO TOP?” 

READ ( 3 , ) M 
C 

10 WRITE! 3, ) "WHICH SP. FREQ. COORD. DO YOU WANT TO FIX, X OR Y?" 
WRITE! 3, ) "TYPE 'X' FOR X, OR 'Y' FOR Y WITH THE QUOT. MARKS!” 
READ ( 3 ) XORY 

WRITE! 5 , ) "SET VALUE FOR THAT SP. COORD. ! " 

READ! 3, 5 FRFIX 
IFCXORt.EQ. 'X' ) 

X=FRFIX 
ELSE 
Y=FRFIX 
END IF 

11 WRITE! 3,) "TO WHAT SP. FREQ. FROM ZERO DO YOU WANT TO PLOT?" 
READ! 3.) FRMAX 

DF=FRMAX/ 100 . 

C 

DO 100 1=1,100 
SPFR! I ) = ! 1-1 ) *DF 
IF ( XORY . EQ . 'X' ) 

Y=SPFR( I) 

ELSE 

X=SPFR! I) 

END IF 

100 CALL FT2D(X, Y,M,TQP , A! 1 , 1 ) ) 

WRITE! 3 , 101 )M, TOP .XORY, XORY, FRFIX 

101 FORMAT!//," N 6 . o£ SEGMENTED TRAFEZOIES , " ,T40 , "M=" , 14 , / , 

+" THE PUPIL FUNCTION IS SPREAD IN" T40,"y = G10.3,/, 

+“ THE " ,A1,"-SPAT. FRE0. IS FIXED AT " ,A1 =",G10.3) 

CALL USPL0 ( SPFR , A, 100 , 100 . 1,1, 13H2D FOURIER TR,13, 

+ lOHS^ATIAL FR , 1 0 , 2HFT , 2 , RANGE , 1H1 , 0 , IE3$ ) 

WRITE! 3,) 

-(-"DO YOU WANT TO PLOT WITH A DIFFERENT SCALE? YES-1. NO-0" 

READ (3.) NY 

IF ( NY.EQ . 1 ) GO TO 11 

WRITE! 3 , 120) XORY , FRFIX 

120 FORMAT!///," YOU HAVE BEEN SETTING " , A1 , " = ", G10 . 3 , " SO FAR.",/, 
+ " DO YOU WANT TO PLOT WITH A DIFFERENT SETTING? YES- 1 , NQ- 0 " , ) 
READ! 3, ) NY 
IF ! NY . EQ . 1 ) GO TO 10 
WRITE! 3?) 

+ " DO YOU WANT TO ANALYZE THE SHAPE WITH A DIFFERENT DPHI?" 
READ! 3, )NY 
IF! NY .EQ . 1 )GO TO 9 
STOP 
END 


C 

r * . 

c" 

c. 

c 


SUBROUTINE FT2D(X,Y,M,T0P,AFT) 

G ! ‘ HUGH ‘ SONG ) ' updated * on APRIL 17, 1986 
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PART II. ERROR ESTIMATION OF SINGLE-COLOR LASER RANGING WITH A 
CUBE-CORNER-RETROREFLECTOR ARRAY AND ITS APPLICATION TO LAGEOS 

ABSTRACT 

The performance of a correlation estimator for the single-color laser 
ranging with a cube-corner retroreflector-array satellite as the ranging target 
is analyzed. The pulse shape of the laser source and the impulse response of 
the photodetector being considered as Gaussian shape, pulse broadening by split 
reflections at the array of cube-corner retroreflectors is analyzed. In this 
case of Gaussian approximations on the source and the detector, an expression 
for the root-mean-square of the timing error in ranging has been obtained under 
the assumption that shot noise and speckle-induced noise are not severe. Such 
analysis has been applied to a simulated ranging experiment with an existing 
satellite LAGEOS having 426 cube-corner retroreflectors to estimate the lower 
bound for the variance of the timing error. 
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I. INTRODUCTION 

An array of cube-corner-retroreflectors (CCRs) has been used as a target 
for laser ranging. It may be employed in land surveying. It has long been 
used in satellite ranging in the National Aeronautics and Space Administration 
(NASA) [1], By measuring the roundtrip time of pulsed laser light, we can 
estimate the distance from the ground station to the target. Among various 
methods for detection and estimation of received light signal, the correlation 
algorithm [2] is studied for ranging in this report. The correlation 
estimation technique has several advantages over other techniques. That 
considers both the shot noise and time-resolved speckle. It does not introduce 
a signal dependent bias that appears in the leading edge threshold technique 
[3]. It is optimum even for non-Gaussian pulses. All current laser ranging 
systems use one wavelength (single-color) of laser light. Hence analysis on 
single-color ranging is more needed than that on two-color ranging. In fact, 
the correlation algorithm for single-color ranging is simpler to analyze than 
that for two-color ranging. 

In this report, explicit formulas are presented for single-color ranging 
with the target having an array of identical CCRs whose scattering cross 
sections may vary according to their locations on the target. The algorithm is 
then applied to a simulated ranging experiment with the Laser Geodynamic 
Satellite (LAGEOS). This satellite plays a key role in NASA's Earth and Ocean 
Dynamics Application Program [1], It was designed as a passive long-lived 
target with a stable well-defined orbit. It has 426 CCRs of a circular-front- 
face type, so that laser light may be reflected and returned to the ground- 
based station. Variation of the scattering cross section according to the 
location on the LAGEOS has been studied. The root-mean-square (RMS) of the 
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ranging time error is calculated analytically and discussed with some exemplary 
situations of the LAGEOS . 

II. THE CORRELATION ALGORITHM FOR SINGLE-COLOR RANGING WITH AN ARRAY OF CUBE 
CORNER RETROREFLECTORS 

The pulsed laser light beam is assumed to have a sufficient beamwidth and 
a homogeneous wavefront over the surface of the target when it reaches the 
target surface. In this incidence we shall consider the signal amplitude at a 
certain point on the target, e.g., the point nearest to the light source, as a 
sequence of pulses. Individual pulses are well separated in time from each 
other and the RMS amplitude of those pulses is normalized and represented by 
f(t). After the light is reflected at each CCR, it is subject to diffraction. 
In addition, the signal power P(t) being detected at the receiver has 
uncertainty since the light components reflected at different OCRs have certain 
path differences between those components and therefore interfere with each 
other unless the individual pulse is short enough to resolve all CCRs 
spatially. We assume that the CCR array changes its aspect with respect to the 
observer slowly, but fast enough to cancel the net interference effect in the 
long run. The power averaged over such random phases can be expressed as 

M 

<P(t)> = A$ ajf(t - Y m )| 2 (1) 

m=l 

without any cross terms. Here A is a proper proportional constant. The 

summation is taken over all CCRs that have nonzero lidar cross sections a . Y 

m m 

is the time dilation associated with the reflection at the mth CCR. Taking an 
average over the random phase, one may express the covariance of P ( t ) as 
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where the asterisk, refers to complex conjugation. The signal-to-noise ratio 
(SNR) is then defined as 


R SN = 


IX <P(t)>dt] 2 


jx Cpdj 

00 


^ 2 ) d t jd 1 2 


(3) 


In the case of direct detection, the mean and the autocovariance of the 
signal at the photon counter output become, according to the Campbell's 
theorem [4], 


<S(t)> = <P(t)> * h(t) ’ 


(4) 


c s (t r t 2 ) = ; P(T) h <V T) h ( t 2" T) dT 


W OD 

+ ( Rw )2 * f C p (T l’ T 2 ) h < t i- T i) h(t 2 -x 2 ) dTj_ dx 2 , (5) 


where h is the efficiency of the photodetector and hw is the energy of a single 
photon. h(t) is the normalized photodetector-impulse-response which may be 
approximated as a Gaussian function, t), with a RMS pulsewidth, a^, 

where 


G(a, t) = — - exp(- I--) 
/2n a 2 a 


( 6 ) 
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If we also approximate the transmitted pulse intensity as another Gaussian 

2 

function, i.e., |f(t)| = G(a^, t), then <S(t)> and C g (t^, t 2 ) in Eqs.(4) and 

(5) become [5] 


<S( t)> = <N> f (t, a ) , 

b O 


(7) 


l l +t 2 


c s (t r t 2 ) = <N> G(/?a h , V t 2 ) f s (— , a s ) 


+ ^7 g( ^V VV f sp<^>’ 


( 8 ) 


where <N> is the average number of detected photons per pulse, and 


r 2 2 
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According to Eq.(3), the SNR is found for a uniformly illuminated array of CCRs 


'SN 


M 

(l 

m=l 
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A correlation estimator is used to determine the arrival time of the 
reflected signal. A schematic diagram is depicted in Fig. 1 to describe the 
principle of the correlation estimator for single-color ranging. A laser pulse 
with the pulsewidth a^ (assumed to be a non-random variable) deflected by a 
beam splitter is detected at the detector directly at the ground station. 

The width and timing of the laser pulse may be measured at this stage. Here 
most of the pulsed laser power is directed to the target (LAGEOS). Being 
reflected at the target, the light returns to the ground station and is 
detected at the detector D 2 . The pulsewidths of both detectors' impulse 
responses are the same, and are a^. Let x be the random variable representing 
the time elapsed until the pulse returns from the target. Then x is to be 
measured from the two outputs of the detectors and D 2 * These output signals 

are supposed to be processed with infinite speed for correlation. 

Each output signal S(t) at D 2 may be quite different from the previous 
ones since the target changes its aspect with respect to the light source even 
though the aspect is very slowly changing. This continuous change of S(t) is 
due to speckle. The speckle-induced noise can change the detected signal 
pattern because individual light pulses reflected at different CCRs interfere 
with each other unless light pulses are narrow enough to resolve individual 
CCRs. 

Now we shall suppose an expected received signal S(t) by simulation as if 


Figure 1. (a) Schematic diagram of the single-color laser ranging system with 
LAGEOS as the ranging target, (b) Pulse patterns at the output of 
the photodetectors D. and D„ with a Gaussian-shaped impulse response 
(FWHP = 2^) and a Gaussian-shaped pulse laser having FWHP = 2a^. 

In the same time scale the expected pattern S(t) is also 1 

illustrated. See text for details. 
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the LAGEOS were laid right behind the beam splitter and illuminated by a 
uniform plane wave. We then correlate S(t) with the actual received signal 
S(t) as 

00 

R(t) = ; S(t-t) S(t) dt . (13) 

CD 

This is the cross correlation function of the two signals when S(t) is assumed 
to be stationary. Here we assume that, compared to the motion of the target, 
the repetition rate of the light pulse source is so high that we may neglect 
the change of the aspect of the target with respect to the observing station. 
Still, it is assumed that the slight change of path-length difference of light 
components reflected at different CCRs gives rise to totally random phase- 
difference between those light components. In this respect, we can use the 
expectation as the average over the random phase. Also we assume that there 
exists an accurate delay to locate the expected signal S(t-T) when correlating. 
Then the correlation estimator can be written mathematically as 


X COR = ar E>l max R(T)] . 


(14) 


Using a Taylor-series expansion [6], we linearize dR(r)/dt around T = T , the 
representative round-trip time obtained by astronomical computation or by 
averaging past several t C q R ' s * 


d 

dt 


R(t) = ^ R(t) 


+ (T-T ) —j R(T) | 

T r dT 1 T 

r r 


(15) 
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dR(x)/dp vanishes at T = x CQR where R(t) 


is maximum. 


Therefore 


T C0R ~ T r 


dR( t) /d 2 R(x) 

lx /dx 2 lx 
r / r 


(16) 


2 2 

provided that d R(x)/dx does not vanish at x = x . This gives the newly 

measured Note that the validity of Eq.(16) strongly depends on the 

validity of Eq.(15). So when the detector output has rugged peaks due to 

split reflections at CCRs, we should obtain by choosing the maximum among 

several peaks of R(x) according to Eq.(14). 

To evaluate the performance of the correlation estimator, we shall 

consider the mean square error of x nriD when x is ideally chosen at x , the 

uuk r o 

true round trip time. Under the assumption that the shot noise and the speckle 
are not so severe that 


var( lVo > « >2 , 

dx z 'x n dx z >t 

0 0 


(17) 


the variance of x^^, <Ax^ r > of the correlation estimator, defined as 


<At C0R > ' <(t C0R <T C0R >) > ” 


(18) 


can be calculated. Now we shall assume that by monitoring the LAGEOS we have 
the complete information on the aspect of the LAGEOS - the information on which 
face directs the ground station with what azimuth angle. Then, since both the 
speckle and the shot noise induced fluctuations in one shot of detected pulse 
is uncorrelated with those of previous pulses, we may assume that supposed 
expected detected signal may have the structural peaks due to timely resolved 
CCRs and that S(t - t ) is equal to <S(t)>. Therefore 
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This enables us to express the bias <t cqr > - t q of the correlation estimator in 
terms of the mean function of S(t) 


<t cqr > - T q = J <S(t)4<S(t)>dt / J <S(t)>^ T <S(t)>dt 
’ ~ d t 


( 20 ) 


As long as <S(t)> is a normalizable function, the denominator vanishes. 

Consequently, <t cqr > * ^ 0 * This result is quite reasonable since is the 

true round trip time. The same reasoning applies to derivation of the 

2 

expression for <At cqr >. 



3 3 

/ _ i ytj.tj) ^-<s( tl) > — <s(t 2 »dt 1 dt 2 


-CD -CO 


( 21 ) 


[J <S(t)> ^-«<S( t ) >d t ] 2 
3t 


As mentioned in the case of two-color ranging [2], Eq.(21) gives the lower 
bound of the variance in the correlation estimator. The actual timing variance 
is expected to be higher depending on the actual SNR and owing to nonideality 
of the delay. Besides, lack of exact information on the aspect of LAGEOS makes 
us to suppose only a smoothed expected pulse pattern S(t), which makes the 
actual timing variance in the case of single-color ranging greater than that 
estimated by Eq.(21). Then if we assume <N> is so large that we can neglect 
the first term of C s ( t ^ , 1 2 ) in Eq.(8), <At^ 0R > is found to be approximately 
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On the other hand, if <N> is not so large, we should retain the first term of 
Eq.(8) into the integral in Eq.(21). Under the assumption that a^ » a^, one 
can reduce the integral due to the first term into the following summation. 
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to (22). Then we should add this to (22) to evaluate > when <N> is not 

so large and cr^ >> cr^. 


III. APPLICATION TO LASER RANGING WITH THE LAGEOS 

Once the information on the lidar cross sections and mutual optical path 
difference of the CCRs is given, one may readily find the root-mean-square of 
the timing error of ranging by summation only. Study on lidar cross sections 
of a CCR array has been done by the present authors [7]. When one estimates 
the lidar cross section of a CCR on the LAGEOS, he has to consider the follow- 
ing factors. First, the front face of a CCR is tilted from the normal 
illumination. Second, the triple-corner-mirrors are not coated.. Third, the 
laser light is linearly polarized. Fourth, there are some manufacturing 
errors including an intentional error in the dihedral angle. However, an exact 
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theoretical estimation of the lidar cross section formula in which all the 
above facts are considered properly is not yet available. As we mentioned in 
Ref. [7], the intentional error in dihedral angle and other manufacturing errors 
make the computation of the lidar cross section very difficult. Numerical 
estimation formulas for such computation was found by the present authors in 
Ref. {7]. Hence empirical evaluation of its variation deserves our notice [8J. 
However, for evaluation and analysis of various detection methods for ranging 
with CCR arrays, a reasonable numerical formula for estimating the lidar cross 
section of a CCR is needed. 

For analysis and evaluation of the timing error, we shall use the results 
in Ref. [7], in which the above-mentioned first three factors for estimating the 
lidar cross section were considered in the computer program completed for the 
first two factors. So in the computer program we are using here, it is assumed 
that the corner-mirror system is not coated, but it is neglected that the six 
sectors of the effective reflecting aperture of one CCR reflect the six light 
components that interfere with each other and consequently make an interfer- 
ence pattern at the receiver site decreasing the lidar cross section irregular- 
ly. So the following results should be revised when we incorporate the exact 
estimation formula in Ref. [7] for the lidar cross section when linearly 
polarized laser light is illuminated into the computer program. 

Here we show in Figures 2, 3 and 4 the expected pulse pattern when we do a 
ranging experiment with LAGEOS. In Fig. 2, the functions f s (t) and f g p(t) have 
been plotted with = 20 psec and with an aspect of the LAGEOS whose 

south pole is directed to the ground station. The pattern depends on the 
aspect of the LAGEOS. Pulse patterns in Fig. 3 were obtained with varying 
pulsewidth with this aspect of LAGEOS. In this case, only 20 CCRs out of 426 
can retroreflect by total reflection at the triple mirror system, while 87 CCRS 



0.20 


The normalized pulse pattern f (t) (solid line) and its speckle 
function pattern f (t) (dashed s line) with a, = cr, = 20 nsec when 
the south pole of lRgEOS directs the ground station. 



0.40 


TIME (nsec) 
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TIME (nsec) 


Figure 3. The expected pulse patterns S(t) or f (t) that would be obtained 
after normalizing to each maximum and averaging over the random 
phase due to speckle and after removing the shot noise of the 
photodetector. Patterns are drawn with varying pulsewidth Or with 
fixed = 0 when the south pole of LAGEOS directs the ground 
stationror the purpose of comparison with those in Ref.[l]. 
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TIME (nsec) 


Figure 4. The expected pulse patterns S(t) or f (t) that are obtained with Q 
varying pulgewidth a, = a, when a LAGloS coordinate (longitude 10 
latitude 85°) directs to the ground station. According to the order 
of degree of ruggedness, the patterns are obtained with = a, = 2, 
5, 10, 20, 40 psec. The most rugged pattern = 2 psec; 

shows the spatially resolved OCRs with a few exceptions at 0.277 and 
0.33 psec. 
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can do so with positive effective retroreflecting areas. Those in Fig. 4 were 
obtained when the aspect of LAGEOS is quite arbitrary (latitude 85° and 
longitude 10° on the LAGEOS coordinate) but still fixed with 

respect to observer. In this case, 23 CCRs can retroreflect the incident light 
by total reflection, while 99 CCRs can do with positive effective 
retroreflecting areas. In computing the lidar cross section of each CCR, the 
velocity aberration effect due to the LAGEOS linear velocity on its orbit has 
already been considered. That is, the maximum spot of the far-field 
diffraction pattern was assumed to lie 34 yurad off-axis from the ground station 
(1]. The former aspect is one of the two aspects of the LAGEOS (that the north 
of the south pole faces the ground station) with which the light components 
reflected at different CCRs interferes most severely. On the contrary, the 
second aspect of the LAGEOS shows a pulse pattern which is virtually 
uninterferred when the pulsewidth of the laser light and that of the impulse 
response of the detector is narrow enough to resolve every pule reflected at 
each CCR. This feature can also be inferred from the timing error analysis. 

The SNR and the RMS timing error for each case hae been calculated 
according to Eqs. (12) and (22) under the assumption that <N> is sufficiently 
large. In Table 1, they are tabulated for the two foregoing aspects of LAGEOS. 
It was found that, except almost improbable particular aspects, like the first 
case of the LAGEOS aspect (south or north pole) and like that of Ref. [2] (a 
flat array of CCRs whose faces direct the laser light exactly), the timing 
error can be greatly reduced up to the quantum limit when we use extremely 
short light pulses and an extremely fast photodetector. The RMS values of the 
timing error in Table 1 were obtained without considering the fundamental limit 

_3 

imposed by wave mechanics. Therefore values less than 10 psec are 
theoretical values when the light source with a wavelength much shorter than 
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TABLE 1. The root-mean-square values of timing error (^ CO n) and the signal- 
to-noise ratio (R_„) for the two aspects of the laGEOS with varying 
(Sc, and o^, the pulsevidths of the laser and the photodetector, 
respectively. 


(a) For 

the aspect that 

the south pole faces 

the ground station. 

<r f (psec) 

<r h (psec) 


AT C 0 R ( pse c) 

r sn 

0.1 

0.1 


5.431xl0" 5 

4.4710 

1 

1 


5.431xl0 -5 

4.4710 

3 

3 


8.420xl0" 5 

4.4710 

10 

10 


6.390xl0 -1 

4.3730 

20 

20 


8.08 

4.020 

62.5 

0 


32.98 

2.0619 

125 

0 


46.07 

1.444 

250 

0 


58.98 

1.172 

500 

0 


64.29 

1.074 

1000 

0 


65.54 

1.046 

(b) For 

the aspect that 

the LAGEOS coordinate 

- latitude 85°, 

longitude 

10 - faces the ground station. 


a f = 

*h 

AT C o R ( Psec) 

r sn 

0.1 


8.299xl0~ 7 

1006.1 
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6. 321xl0 -2 

83.263 
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0.1283 


59.382 

5 


2.1235 


20.856 

10 


6.2758 


9.4821 

20 


17.040 


4.7458 


40 


25.161 


2.6683 
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the visible range of wavelength is used. In Appendix A, a computer program is 

given, which computes the SNR and the RMS timing error and plots the expected 

pulse pattern f (t) and the speckle pattern f (t) for a given aspect of LAGEOS 
s s p 

and given pulsewidths of the laser and of the photodetector. 

CONCLUSIONS 

In this report, the correlation estimator for laser ranging with a 
single-color laser has been studied. Pulse patterns which will be obtained by 
actual ranging experiment are obtained by simulation with the RMS timing error 
that will accompany in somewhat idealized experiment. It is idealized because 
we usually do not know the aspect of the CCR array unless it is fixed. 

Gaussian approximation has been employed throughout the analysis. That is, the 
laser pulse and the impulse response of the detector have Gaussian pulse 
shapes. An expression for the variance of timing error with a given array of 
CCRs characterized by the lidar cross section and the mutual pathlength 
differences was found. It was then applied to LAGEOS with somewhat incomplete 
data on lidar cross sections of each CCR. It is incomplete since the 
polorization change has not been simulated. It should be revised when the data 
for correct lidar cross sections is available. 

Aside from incomplete datas for the lidar cross section, the timing error 
analysis has been based on the fact that the expectation of the received signal 
is fully known. For this statement to be true, we must have complete 
information of the aspect of LAGEOS, which may hardly be obtained. Therefore 
the timing error analysis made in this report gives the theoretical limit or 
lower limit of the timing error. However this work will provide the basis for 
comparison with other detection techniques and will also provide the way of 
evaluating various CCR arrays that were and will be developed for ranging 
application. 

The support from NASA for this work is greatly appreciated. 
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APPENDIX 


>P 

PROGRAM VEL 


G. HUGH SONG updated on FEB. 21, 1987 
Plots the signal pulse shape function Fs(t) 
and the speckle pulse shape function Fsp(t).. 
along with the root mean square error of ranging. 

The velocity aheration is considered in this program. 
Finally, the rms ranging error (RMSERR) is calculated. 

ASPECT ANGLES : INCIDENCE AND POLARIZATION ANGLES OF LIGHT 
IN THE LAGEOS COORDINATE SYSTEM 
ALA (LATTITUDE), ALO < LONGITUDE) , DPCN IN. dag. 

H ( HORRIZONTAL ) , V (VERTICAL), FON IN radian. 

THE VELOCITY ANG LE 

DVDN (deg), VDNtrad) : MEASURED W.R.T. THE NORTH POLE 

THE LIDAR SYSTEM PARAMETERS ARE DEFINED IN STATEMENT #302 
PHI, GAM, & POL ARE DEFINED IN SUBROUTINE SUBCCR 

DIMENSION RANGEC4) ,FSY< 201,1 ),F3PY( 201,1) 

COMMON /PULBLOCK/sfGFdO) .SlGn(lO) ,SNR(10) ,RMSERR(10) , 

+ NOCCR , NEFCCR , PTHMI , FTHMA , 

+ R , C , COVTs IV , H , COH , S IH , LFOL , LFRESN , PON . VEL . VDN 

COMMON /CRABLOCK/T( 201) ,F3( 10,201) ,FSP( 10,201) 

DATA C/2.9979E3/ 

RAD ( DEG ) = ATAN ( 1 . )*DEG/45. 

CALL LAGEOS (R) _ _ _ ‘ 

WRITeTSTT" INPUT THE PLOTTING TIME RANGE TI £ TF IN nsec." 
READ ( 3 , )TI ,TF 
TI=TI*l.E-9 
TF=TF*l.E-9 


NG=0 

NG-NG+1 


[GO TO loop count 


IF ( MG . EQ . 1 ) GO TO 3 

WRITE (3,)" WILL YOU USE THE SAME TIME SCALE? Y-l, N-0!" 

READ (3.) NY 
IF(NY.EO.O) GO TO 1 

FORMAT? “ ~Will you use the same ASPECT OF THE LAGEOS,",/, 

+ •■ the same VELOCITY DIRECTION ,",/ , „ 

+" and the same POLARIZATION OF THE LASER? Yes-1, No-0") 

READ( 3 , )NY 

IF(MY.EQ.l) GO TO 204 

FORMAT? / / " 1 INPUT THE ASPECT ANGLE OF THE LAGEOS'",//, 

+" THE LONGITUDE (DEGREES) (0. -< ALA =< 360.)") 

PEAD( 3 ) ALA 

WRITE ( 3 , ) " THE LATTITUDE (DEGREES) (-90. =< ALO < 90.)" 

READ( 3 . ) ALO 
H=RAD( ALA) 

V =RAD ( ALO ; 

COV=COS(V) 

5IV=SIN( V) 

COH=C0S ( H ) 

3IH=3IN ( H ) 


WRITE ( 3, ) 

+" INPUT THE VELOCITY DIRECTION W.R.' 
READ ( 3 , ) DVDN 
VDN =RAD ( DVDN ) 


THE NORTH POLE IN deg ! " 


FORMAT? 3 ARE YOU GOING TO CONSIDER: THE FRESNEL REFLECTION",/, 
+" INSIDE THE CCR AS WELL AS THE TOTAL REFLECTION?",/ 

+ " It will take 20 mir./set of plots. Yes - a, No-u . ,/) 

READ (3,) LFRESN 
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FORMAT ( 1 , 1 ." IS THE INCIDENT LIGHT POLARIZED? YES-1, NO-O") 

READ (3,) LPOL 

IF ( LFOL >203,204,203 

WRITE ( 3,308) 

FORMAT!/" Sorry. Not developed for polarized light."/. 

+ " The procrram assumes unpolarized liarht in the f oliowmcr. " ) 

j." r.Ttr atsr.r r acciTur 'Tn_-7i'T> T-iic r rr'i.rr to rrunrnorv Dnr acr7m n , 


+,/," DIRECTION OF THE LAGEOS.") 

READ(3,)DP0N 
PON =RAD ( DPON ) 

WRITE! 3 ,302) 

FORMAT (//," INPUT THE LIDAR SYSTEM PARAMETERS!",// 

+" SIGF(sec) : RMS LASER PULSE WIDTH,",/, 

+ " SIGH(sec) : RMS WIDTH OF THE RECEIVER IMPULSE RESPONSE.",/// 
+ " How many diff. set of SIGF Si SIGH do you want? JC = ?") 

READ (3,) JC 

WRITE (3, 309) 2.E9*R/C 

FORMAT!//" Note that 2*P.adius/C = " .G13.6," nsec.",/, 

+ " Input SIGF and SIGH in psec line by line for each set!") 

DO 229 J=1,JC 

READ! 3 , ) SIGF! J) , SIGH(J) 

SIGFfu )=SIGF( J)*1.E-12 
SIGH! J)»SIGH( J)*1.E-12 

DO 55 IT-1,201 

T! IT) =TI+( IT-1 ) / 200 . * (TF-TI ) 


-DO LOOP FOR SETS OF Fs(t) and Fsp(t)-- 


DO 9993 J»1,JC 
CALL PULSEFAT ( J ) 
CONTINUE 


of the DO LOOP 


WRITE! 77 , ) JC ,NCCCR .NEFCCR , PTHMI ,PTHMA 
DO 175 J=1,JC 

WRITE (77,) SIGF! J) ,SIGH(J) , SNR ( J ) , RMSERR ( J ) 

RANGE ( 1 5 =T ( 1 ) 

RANGE ( 2 ) =T ( 20 1 ) 

RANGE! 3) =0. 

RANGE ( 4 ) = 0 . 

WRITE(3,)" J = what d'yu plot on screen? If J=0, no plot." 
READ(3,) J 

WRITE! 5,352)1. E12*SIGF( J ) , 1 . E12*SIGH( J ) 

FORMAT! IX, "LASER PULSEWIDTH : S IGF= " G16 . 7 "nsec . " , / , 

-" RECEIVE?. IMP. RES? . PULSEWIDTH : SIGn= " ,Gl£ . 7 , "psec . A , / / ) 
IF(J.EQ.O) GO TO 120 

DO 119 1=1,201 
FSY! I,1)=FS! J,I) 

F3FY (1,1) =FSP ( J , I ) 

CALL USPL0 ( T , FSY ,201,201,1.1. 5HFs ( t ) , 

4 , 4HTIME , 4 , 5HFs ! t ) , 5 ,RANG±,1H1 , 0 , IER) 

CALL USPLO ( T , FSPY ,201.201,1, 1 , 6HFsp ( t ) , 6 , 

+ 4HTIME , 4 , SHFsp ( t ) , 6 , RANG?: , IHi; 0 , IER ) 

WRITE (3,51) RMSERR ( J)*1.E12, SNR ( J ) , NEFCCR , NGCCR , 1 . E9* PTHMI , 1 . E3 
+ *FTHMA 


■" <-- not valid when you done consider Fresnel ref.",/ 
■" from PTHMI = " ,G10 . 3 , nsec to ?THMA=" ,G10 . 3 , "nsec . " , / / 
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C 

c 


WRITE(3,)" Another plot in the screen? Yes-1, No-0.” 

READ t (3, ) NY 

379 wrjtE( 3? )” ^ Destroy current dath knd get new? Yes#l, No-0." 
READ( 3 . ) NY 
IF(NY.GE.l) GO TO 2 

WRITE ( 3 , ) " WANT TO PLOT ON THE ELECTOSTATIC PLOTTER? Y-l,N-0" 
GO TO 999 

31 FOR^Tf r ' 3 WHAT DO YOU WANT TO PLOT?",/ w 

+" 1 set of Fs ( t ) and Fso(t)? Then type No. o^ the set n J! J , 
+ Several ( .LE. ten) Fs(tT s on the same sca^e. Type 11. 

+ " OR do you want to stop? Type 99 J l ) 

SftS.LE.'lo') CALL GRAPHICS (J) 

IF(J.EO.ll) CALL MULTIGRA(uC) 

IF( J.EQ. 999) STOP 
GO TO 379 
399 STOP 
END 


C 


c 

c 

c 

c 

c 

c 

c 

c- 


SUEROUTINE PULSE? AT ( J ) 


XSC ( m ) 
3X3 C 
SPG AH 
SGAM 


CROSS SECTION OF THE m-TH CCR 
TOTAL SUM OF X3C(m)'S 
I N T ERFER ENCE FACTOR FOR EAch CCR 
WEIGHTED SUM OF SPGAM'S 


;io) , 


DIMENSION XSC ( 426) , PTH ( 426) .... 

COMMON / GRA3L0CK / T ( 2 0 1 > ,FS( 10,201) ,FSF(10,2ul) 
COMMON / PULBLOCK / SIGF ( 1 0 ) , S IGH ( 1 0 ) , SNR ( 1 0 ) , RMSERR ( 

COMMON / BOARDER / PHI , TAP 
DATA XSC, PTH/ 35 2*0./ . 

DATA Disl 1 , VP / 5 . 9E+b , 34.E-6/ 


C 

C 

c- 


DATA RlTVT^DiAlWL/ I ’ 455 , .0273, .0331, 5.32E-7/ 

DATA NOCCR .NEFCCR/ 2*0 / 

RAD ( DEG ) =ATAN ( 1 . )*DEG/45. 

?I=4.*ATAN( 1. ) 

SQSG4=4 . * ( SIGF ( J ) **2+SIGH( J ) **2 ) 

SlGG=3QP.T( SQSG4 ) / 2 . _ ... 

Compute Fs ( t ) . 

PTHMI = R / C 
PTHMA=0 . 

3XSC=0 . 

3GAM=0. 

DO 81 N=1 , 426 
RTHE=RAD < THE ( N ) ) 

RPSI=RAD( ?SI ( N) ) 

RORI=RAD( ORI ( N) ) 

COPHI=vCOV*(COH*X(N)+SIH*Y(N) )+SIV*Z(N) ) /R 
IF ( C0PHI . LE. 0 . ) GO TO 30 
IF ( COFHI . GE . 1 . ) PHI = 0 . 

IF ( C0FHI . LT. 1 . )FHI=AC0S( COFHI) 

GAM=ATAN(CGV*SIN(RPSI-H> / ( SIV*SIN i P.THE) -C0V 
+ *CQS ( RTHE) *CCS ( RPSI-H) ) ) -RORI 

2ET = ATAN ( 3 IN ( P.THE ) * 3 IN ( h’-RPS I ) / ( COS ( P.THE ) *L 0 V- S IN t RThr. ) 
+ ^ *SIV*C0S ( H-RPSI ) ) ) 

?0L=PCN-Pl / 2 ! +ZET ! TE ANGLE FOR EACH CCR (POL=0 : TE) 

F0L=PI/4 . 
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END IF 
VEL=VDN-ZET 

C ! VEL : VELOCITY ANGLE FOR EACH CCR(VEL=0 : SHORT AXIS=x-AXIS) 

Sc-|k*^S*? t 5eK‘ /DIST/DIA !VP : AEERATION ANGLE! rad) 

YK=RK*SIN ( VEL ) 

CALL SUBCCR(1.455.GAM,FHI,0,0,LFRESN,POL,ALP,BET, 

+ DFHIM,ARN,TOTREF,REF,TOP) 

c 

IF(REF.EQ.O. ) GO TO 80 
NOCCR=NOCCR-t-l 

IF ( TCTREF . GT . 0 . ) NEFCCR=NEFCCR+1 
C \ Count the No. of effective CCRs. 

CALL FT2D(XK, YK, . 01 , TOP, FT) 

XSC(N)=FT*FT*REF ! Velocity aheration effect 
PTH (N) =2.*(R*( 1. -COFHI ) +VT*SORT ( RI *RI - 1 . +C0FHI *COPHI ) ) / C 
C i REFLECTION TIME CORRECTION FOR EACH CCR 

FTHMI = AMIN1 ( PTHMI , PTH ( N ) ) 

FTHMA* AMAX1 ( PTHMA , PTH (N) ) 

C \The first and last CCRs' locations in time. 

DO 98 IT=1 , 201 
TIPN=T( IT) -PTH(N) 

93 FS ( J , IT ) *FS ( J , IT ) +XSC ( N ) *EXP ( -TI?N*TIFN / S IGG* * 2 ) / S IGG / SQRT ( PI ) 

GO TO 81 
8C XSC (N) =0 . 

31 SXSC=3X5C+XSC(N) 

C 

C Compute Fsp(t) the the timing error. 

RMSDEN* -XSC ( 1 ) *XSC ( 1 ) / 2. 

RMSN0M=0. 

NACC=1 

C 

DO 85 N=2 ,426 

IFCXSC(N) .EQ.O. )G0 TO 85 

NACC=NACC+1 

WRITE! 3 . ) J.NACC.N 

RMSDEN*RMSDEN-XSC ( N ) *X3C ( N ) / 2 . 

NM1=N-1 

DO 84 M*1 ,NM1 

IF ( XSC ( M) . EQ . 0 . ) GO TO 34 

TD=FTH(M) -PTH(N) 

TM»(PTH(M)+PTH(N) ) / 2. 

SPGAM=XSC(M)*XSC(N)*EX?(-TD*TD/4. /SIGFf J)/SIGF( J) ) 

3CAM=3GAM+2 . * SPG AM 

RMSDEN =RMSDEN+2 . *XSC ( M) *XSC (N) * ( TD*TD/ SQSG4- . 5 ) 

+ *EXF(-TD*TD/SQSG4) 

DO 37 IT=1 , 201 
TITM*T( IT) -TM 

87 FSP ( J , IT) =FSF( J , IT) +2 . *SPGAM 

+ *EXP( -TITM*TITM/3IGG/SIGG) /SIGG/SQRT! PI ) 

DO 88 NP* 1,426 
IF(XSC(NF) .EQ.O. )G0 TO 88 
TFN*TM-PTH!NP) 

CSNOM=0 . 

C 

DO 89 MP=1,NP 
IFCXSC(MP) .EQ.O. )G0 TO 89 
TFM=TM-PTH ( MP ) 

SNOM=XSC ( MP ) *T?M*EXP ( -TPM*TPM/ SQSG4 ) 

IF ( NP.EO . MP ) SNOM=5NOM/2 . 

CSNOM*CSNOM+SNOM 
89 CONTINUE 

RMSN0M=RMSNQM+CSNQM*2 . *XSC ( NP ) *TFN*SFGAM 
+ *EX? ( -TFN*TPN / S QSG4 ) 

33 CONTINUE 

34 CONTINUE 
55 CONTINUE 


SNR ( J ) =SXSC*SXSC/SGAM 
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339 


97 

C 

83 

C310 

C 


c 

c= 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 


RMSN0M=RM3NQM/ SGAM/ 4 . 

RMSDEN = RMSDEN / SXSC . 

RMS ERR ( J ) =ABS ( SORT ( RMSNOM/ SNR ( J ) ) / RMSDEN ) 

Normalise Fs(t) and. Fsp(t). 

DO 97 IT-1,201 
FS ( J , IT ) =FS ( J , IT ) /SXSC 
FSP ( J , IT ) =FS? ( J , IT ) / SGAM 
DO 83 1*1.201 

WRITE ( 66 , 310 )T( I ) ,FS ( J ,1 ) Tx 

WRITE ( 66 f )T( I ) ,FS( J , I) ,FSP< J,I) 

FORMAT ( 2 f IX , G1 6 . 7 ) ) 

RETURN 

END 


SUBROUTINE LAGEOS(R) 

PRODUCES SIX ARRAYS OF COORDINATES AND ORIENTATIONS 
AND THE COMMON RADIUS OF ALL 426 CCR'S. 

X.Y.Z: CARTESIAN COORDINATES 
SPHERICAL COORDINATES 

R : THE COMMON RADIUS _ _ 

THE: LONGITUDE (NORTH POLE : THE=0 . ) 

PSI: LATITUDE 

ORI : ORIENTATION OF ONE EDGE OF A CCR 


DIMENSION DTHEC20) .DORI( 20) ,BETW( 20) , NUM 20) 
COMMON /LAGE/X( 426 f ,Y(426) ,Z ( 426 ) ,THE( 426) ,PS*^~w 
DATA NUM/ 1,2,3,20,33, 61,83.119.150,182. 

+ 2l4',246, 178:309:346, 367,390,468,420,426/ 


PSI ( 426 ) , ORI ( 426 ) 


DATA 



94 

140.691, 


150.421, 160.151, 169.381, ISO./ 


97. , 63. , 55. , 47. , 
78., 90., 18., 102./ 
7*0 . / 


341: 
342: 
343 : 
244 : 
345: 
346 : 
347: 
343: 
349: 
350: 
351: 
352 : 
352 : 
354: 
355 : 
356: 
357: 
358 : 
359: 

360 : 

361 : 
262: 
363 : 
3 64 : 
2' 65: 

o 66 : 

367 : 
368: 
369 : 


383 

C 

339 

10 

20 

387 


C 

C = = = 
C 


DATA DORI/55. , 91. , 43. , 31., 29., 11.. 

+ 94. ,102. ,110. , 24.. 58., 76., 

DATA SETW/ 7*0. , 1., 0., 1., 1., 0., 1., 

FI =4. *ATAN ( 1 . ) 

R=. 2721 

DO 887 J-1,426 
DO 888 1-1.19 
NU=NUM( 1*1 ) 

IF ( J . LT . NU ) GO TO 889 
CONTINUE 

THE ( J ) =DTHE ( I ) 

PSI ( J ) = ( J-NUM( I)+.5*BETW( I) )*360. / (NUM( 1+1) -NUM(I) ) 
ORK J) -DORK I) -26. *( J-NUM( I) ) 

IF ( ORI ( J ) . GE . 0 . ) GO TO 20 
ORK J) -ORK JJ+120. 

GO TO XO 

RPsi=PSI( J)*PI/180. 

RTHE-THE ( J ) *PI / 180 . 

X< J)=R*SIN(RTHE)*COS(RPSI) 

Y( J ) =R*SIN(RTHE) *SIN( RPSI ) 

S( J)=R*COS(RTHE) 

RETURN 

END 


C. 

c 


SUBROUTINE FT2D(X, Y,B,TOP,FT) 


C 

c 

c- 


G. HUGH SONG, updated on APRIL 15, 1986 
SLO : SLOPE OF EACH LINEARIZED TRAPEZOID 
TOP : THE TOP POINT OF THE PUPIL 
B : THE HEIGHT OF THE TRAPEZOID 

BP : THE HEIGHT OF THE ISOSCELES 
X, Y : THE SPATIAL FREQUENCY 
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c 

c 


N=IFIX(TOP/B) 

BP=T0P-N*B 
SL0=-F(N*B) /BP 
TMB=T0P-BP 
TMB2=T0F-3P/2 . 

IF(X.EQ.O. ) 

FT= 2 . *SL0*BP* ( TMB*S INK ( 2 . *Y*TMB ) 

+ -TMB2*SINK(2.*Y*TMB2)*SINX(BP*Y) ) 

ELSE 

YPX=Y+X*SLO 

YMX=Y-X*SLQ 

FT=2 . /X*3P*( SIN(TMB2*Y-SL0*BP*X/2 . ) *SINK(BP*YPX> 

+ -SIN(TMB2*Y+SL0*BF*X/2. )*SINK(BP*YMX) ) 

END IF 

DO 99 I-l.N 
BM=(I-.5)*B 

5LG = (F( I*B) -F( ( 1-1 )*B) ) /B 
FM=(F( I*B)+F( (I-1)*B)+2.*F(BM) )/4. 

IFfX.EQ. 0. ) 

FT=FT+4 . *B* ( FM*SINK ( B*Y ) *C0S ( Y*BM ) 

+ +SL0*BM*SINK(2.*Y*BM)*(C0S(B*Y/2. )-SINX(B*Y> ) ) 

ELSE 

YMX=Y-X*SL0 

YPX=Y+X*SLO 

FT=FT+2 . /X*B* ( SIN(X*FM-s-Y*BM) *SINK( B*YPX) 

+ +5IN(X*FM-Y*BM)*SINK(3*YMX) ) 

END IF 

.99 CONTINUE 

RETURN 
END 

FUNCTION SINK(Z) 

IFCZ.EQ.0. ) 

3INK=1 . 

ELSE 

SINK=SIN(Z/2. )*2./Z 
END IF 
RETURN 
END 


C«” = = - 


8= 

c 


FUNCTION F(Y> 

COMMON /BOARDER /PHI, TAP 

F=COS ( PHI ) * ( SQRT( 1 . - Y*Y ) -SQRT ( 2 . ) *TAP ) 

RETURN 

END 


SUBROUTINE GRAPHICS (J) 

DIMENSION XI ( 201) , Y1 ( 201 ) , Y2 ( 201 ) 

COMMON /GRABLOCK/ T( 201 ) ,FS ( 10 , 201 ) ,FSP<10 , 201 ) 
DATA Y2M/0./ 

DO 1 1=1,201 
XI ( I) =T( I)*l.E+9 
Y2 ( I ) =FSP ( J , I ) / 1 . E+9 
Y2M= AMAX1 ( Y2M , Y2 ( I ) ) 

Y1 ( I ) =FS ( J , I ) / 1 . E+9 
WRJTE(_3 f 301) Y2M 


1 

301 


FORMAT ( ' Max of graph=" ,G10 . 3 , " Set it 
+, /." Make it a little larger than the 
£eAD(3,) Y2M 


again for a 
original max. 


nice plot. 
Ymax=? " ) 



Y2M, 6. ,4. ) 


CALL WHERE 
RETURN 
END 


14,1) 
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SUBROUTINE MULTIGRA(NG) 

COMMON /GRABLOCK/X( 201) ,Y( 10,201) ,2(10,201) 

DIMENSION YG ( 2 0 1 ) , YM(5) 

DATA YM/ 5*0 . / 

SCALING 

DO 3 J = I , NG 
DO 3 1=1.201 

YM( J) =AMAX1 ( YM( J) ,Y(J,I) ) 

DO 5 1=1,201 
DO 4 J = 1 . NG 
Y(J.I)=Y(J,I)/YM(J) 

X(lJ«X(I)*I.E9 

CALL PINIT( 1.5, 5., 1000) 

DO 7 1=1^201 

CAl£ ) MGRPAl!x,YG, 201,X(1) ,X(201) , 0., 1., 6., 4.) 
IF(NG.EQ.l) GO TO 999 

DO 6 J=2,NG 
DO 8 1=1,201 
YG( I) =Y( J.I) 

CALL DATPf,2(X,YG, 201. 0,1, 0,0. .0,2,0) 

CALL TITLE ( 1HL , 5HF s ( t ) , 5 , . 14 . 1 ) , . , t 

CALL TITLE ( 1HB , 1 1HTIME ( nsec ) , 11 , . 14 , 1 ) 

CALL WHERE 

RETURN 

END 

SUBROUTINE SUBCCR ( RI , GAM , PHI , LFACE , LCOAT , LFRESN , POL , 

+ ALP, BET, DPHIM , ARN , REF , REFFR , TOP ) 

G. HUGH SONG, FE3. 21, 1987 

COMPUTES THE REFLECTANCE OF A CUBE CORNER RETROREFLECTOR 
FROM SPECIFIED DATA OF THE RI , GAM, PHI ANGLES. 

RI: Refractive index of the CCR. 

GAM ( rad) : Angle of the plane of incidence 

Tfc is measured from one edce direction (north pole) . 



ace 


BET “ANGLE OF' 'PASSIVE ROTATION WITH RESPECT TO THE 2 AXIS 
ALP “ANGLE OF PASSIVE ROTATION WITH RESPECT TO THE Y-PRIME AXIS 
DPHIM=MAX ACCEP TANC E ANGLE OF INCIDENCE TO REACH THE VERTEX 
ARN “NORMAL I ZED EFFECTIVE REFLECTING AREA 
REF=REFLECTANCE FOR INTENSITY 

Be sure to use the P option in SAUF77 compilation^ 

'common /AL3E/SIA,CCP,RT2,SIP,RT3 ,C0A,TAA,?I 


c- 


REF=0 . 

REFFR = 0 . 

ARN=0 . 

RT2=S0RT( 2. ) 

RT3=SQRT( 3 . ) 

PI =4. *ATAN ( 1 . ) 

SIP=SIN( PHI ) /RI 
FH=ASIN ( SI? ) 

!FH=angle of transmitted licrht through the fronc race 
C0P=C0S ( PH ) 

TAP=3IF/C0F 
COCP=CQS ( FHI ) 


514: 

IF ( LFACE. E0 . 

515 : 

TAFM=i. /RT2 

5^5: 

ELSE 

J i / I 

TA?M=1 . /RT2 


/COS ( AMIN1 ( ABS ( PI/ 3 . -GAM) , ABS ( FI-GAM) , AES ( 5 . *PI / 3 , 
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+ -GAM) ) ) 

END IF 
C 

3 IPM=S IN ( ATAN ( TAFM ) ) *RI 
DFHIM=90 . 

IF( SIPM. LT. 1 . )DPHIM=ASIN( SIFM) *180 . /PI 

IF ( TAP. GE.TAPM) RETURN !REF=0. is to be returned. 

C 

CALL ALPBET ( GAM , ALP , BET ) 

C 

C— - — - - 

C COMPUTATION OF AREA 

C 

GAMS 3 AMIN1 ( ABS ( GAM ) , ABS ( GAM-2. *PI/3. ) , ABS ( GAM-4 . *PI / 3 . ) , 

+ ABS ( 2 . *Pl-GAM) ) 

CALL ALPBET (GAMS, ALPS, BETS) 

C 

SIB=SIN( BETS) 

C0B=C0S ( BETS ) 

C 

IF ( LFACE . EQ . 0 ) 

ARN=2. /PI*C0CP* ( ACOS ( RT2*TAP ) -RT2*TAP*3QRT( 1.-2. *TAP*TAP ) ) 
TOP 3 SORT (1.-2. *TAF*TAP ) 

ELSE IF ( LFACE. EO . 1 ) 

AREA=4 . *SIB*COB*CQA/TAA 

IFfTAA. LE. ( 2 . *C0E ) ) AREA=SIB* ( 4 . *C0A-5IA/CC3 ) 

IF( TAA. LE. COB ) AREA 3 SIB*COA* ( 4 . -COB/TAA) 

IF( TAA.LE. SIB) WRITE ( 3 , ) "ERR?? " 

C AP.EA-SIA*( 4. -COB/SIB) 

C IF( (2.*SIB) .LT.COB) AREA=4 . *SIA*SIB/COB 

AP.EA 3 AREA*COCP/COP 
ARN * AREA/ RT3 
END IF 
C 

IF(ARN.LE. 0 . ) 

WRITE! 3,) M Negative ARN??? Error!! ARN = " ,ARN 
ARN*0 . 

RETURN 
END IF 

IFCLC0AT.EQ.1) GO TO 39 

C 

C CHECK THE TOTAL REFLECTION 

COSTH=SQRT( 1 . -1 . /RI/RI ) 

C0MIN=AMAX1 ( COB*COA , S IB*COA, S I A ) 

Pi ( C0MIN * GE . COSTH ) . AND . ( LFRESN . EQ . 0 ) ) RETURN 

C CALCULATION OF REFLECTANCE (REF) 

C NORMALIZED TO THE CASE OF NORMAL INCIDENCE 

39“* TITE= 2 ! *COCP / ( COCP+RI *COP ) 

TITM=2 . *COCF/ ( RI*COCP+COP ) 

T0TE=2 . *RI*COP/ ( RI*COP+COCP ) 

T0TM*2 . *RI*COP/ ( COP+RI*COCP ) 

TMC=SIN(POL) !TM-component 
TEC=COS(POL) ! TE-COMPONENT 
REF 3 ( ( TEC*TITE) **2+ ( TMC*TITM) **2 ) 

C TEC & TMC should be redefined for exit wave in the future 
REF 3 REF* ( ( TEC *TOTE ) * * 2 r ( TMC *TOTM ) * * 2 ) 

REFFR= REF 

IF( LCOAT.EQ . 1 ) GO TO 31 

c 

C and consider Fresnel reflection in triple corner mirror. 

r 

C 

IF( COMIN.GE.COSTH) THEN 
THI =ACOS ( COMIN ) 

C \ Incident ang. to transm. corner mirror 

THT=AS IN ( S IN ( THI ) *RI ) 

C \ Transm. ancr. out of the corner mirror 

RTE=TAN ( THI -THT ) /TAN l THI+THT ) 

RTM=-SIN(THI-THT) / 5 IN ( THI +THT ) 

C \ Fo r p olarized input, TEC and TMC should be redefined in futur 
REFER 3 REF* ( ( RTE*TEC ) * * 2 r ( RTM*TMC ) * * 2 ) 


a^o^cna»cno*iriOio^aka»0'*a)Crt0^a»cncncno'»aiuiuiuiuiuiL , «i».n 
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592: 


31 


c— - 
c 

333 


c 

c== 

c 


c 

c 

c 


600 

700 


C 


REF=0 . 

END IF 
CONTINUE 

jjO NORMAL I ZAT ION w.r.t. 

REF=REF/ ( 4 . *RI/ ( RI+1 . ) / (RI+1. ) )**2 

RETURN 

END 


NORMAL INCIDENCE. 


SUBROUTINE ALFBET ( GAM . ALP , BET ) 

COMMON /ALBE/SIA,C0P,RT2,SIF,RT3,C0A,TAA,PI 


COMPUTE ALP Sl BET l 


C0G=CQS ( GAM) 

SIA=(C0P+RT2*SIP*CCG) /RT3 
ALP=ASIN( SIA) 

C0A=C0S ( ALP ) 

TAA=SIA/CQA 

SB0R2=(RT3*C0P-SIA) /C0A/RT2 
IF ( SBCR2 . GE. 1 . ) GO TO 600 
BET=PI / 4 . +AC0S ( SB0R2 ) 

IF ( GAM. LT. PI ) BET=FI/ 2 . -BET 
GO TO 700 
BET=PI / 4 . 

RETURN 

END 


S IN ( BET ) +C0S ( BET ) =RT2*SB0R2 
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PROGRAM POSTVEL 


G. HUGH SONG updated on FEB. 21, 1987 
Plots the sianal pulse shape function Fs(t) 
and the specJcle pulse shape function Fsp(t). 
along with the root mean square error of ranging. 

The velocity aberation is considered in this program. 
Finally, the rms ranging error (RMSERR) is calculated. 

THIS PROGRAM RESUMES TGHE WORK VEL AFTER A WHILE 


DIMENSION RANGE ( 4 ) ,FSY( 201 , 1 ) .FSPY( 201 , 1 ) 

DIMENSION SIGF( 10) ,SIGH( 10 ) ,SNR( 10 ) ,RM§ERR( 10 ) 

COMMON /GRA3L0CK/T( 201 ) ,FS( 10 ,201) ,FSP< 10,201 ) 

DATA C/2.9979E3/ 

RAD (DEG ) =ATAN ( 1 . ) *DEG/45 . 

READ ( 7 7 , ) JC , NOCCR , NEFCCR , PTHMI , PTHMA 
DO 175 J=1,JC 
DO 174 1=1,201 

READ (66, )T( I) ,F3( J.I) ,FSP( J.I) 

READ( 77, )SIGF( J) ,SIGH( J) ,SNR( J) , RMSERR ( J ) 

RANGE ( 1 ) =T ( 1 ) 

RANGE ( 2 ) =T ( 2 0 1 ) 

RANGE ( 3 ) =0 . 

RANGE( 4 ) =0 . 

WRITE ( 3 , ) " J = what d'yu plot on screen? If J=0, no plot." 

READ ( 3 , ) J 

WRITEf 3,352) l.E12*SIGF( J) ,1.E12*SIGH( J) 

FORMAT ( IX. "LASER PULSEWIDTH : SIGF= " ,G16 . 7 , "Dsec . " , / , 

+" RECEIVER IMP . RESP . FULSEWIDTH :SIGA=" ,Gl6.7, "psec. * , // ) 
IF(J.EQ.O) GO TO 120 

DO 119 1=1,201 
FSY (1,1) =FS ( J , I ) 

FSPY (1,1) =FSP ( J , I ) 

CALL USPLO ( T , FSY ,201,201,1,1. SHFs ( t ) , 

+ 4 , 4HTIME , 4 , 5HF s ( t ) , 5 , RANGE , 1H1 , 0 , IER ) 

CALL USPLO (T, FSPY, 201, 201. 1. 1.6HFsp( t > ,6, 

+ 4HTIME , 4 , 6HFsp ( t ' , 6 , RANGE , 1H1? 0 , IER ) 

WRITE (3.51) RMSERR ( J)*1.E12, SNR ( J ) , NEFCCR , NOCCR , 1 . E9 *PTHMI , 1 . E9 
+ « PTHMA 

FORMAT ( / / / " THE RMS ERROR OF RANGING IS " ,G10 . 3 , "psec . " , / / , 

+" THE SIGNAL -TO -NOISE RATIO IS ",G10.3,//, 

+ " The total No. of effectively reflecting CCRs is ” ,13,//, 

+" The total No. of CCRs with positive ref. area is ",I3,//, 

+" from PTHMI = " ,G10 . 3 , "nsec to ?THMA=" ,G10 . 3 , "nsec .",/// ) 


WRITE(3,)" Another olot in the screen? Yes-1, No-0." 

READ (3,5 NY 
IF(NY.EQ.l) GO TO 377 

WRITE( 3 , ) " WANT TO PLOT ON THE ELECT0STATIC PLOTTER? Y-l,N-0" 
READ (3.) NY 
IF(NY.EQ.O) GO TO 999 

WRITE ( 3.31) 

FORMAT ( r ' WHAT DO YOU WANT TO PLOT?"./ 

+" I set of Fs(t) and Fsp(t)? Then cvoe No. of the set J !",/, 
+ " Several! .LE. ten) FsftVs on the same scale? Type 11!",/, 
t" OR do you want to atop? Type 9991") 

READ (3,! J 

IF(J.LE.IO) CALL GRAPHICS ( J ) 

IF( J.E0.11) 

WRITE 13 , 5 " How many (JC) curves, J=1,...,JC, JC=?" 

READ (3,5 JC 
CALL MULTIGRA ( JC 5 
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/ ** : 
7 d : 
7 6 : 
77: 
73: 
73 : 
30 : 
81: 
32: 
33: 
34: 
35: 
36: 
87: 
33: 
33: 
30 : 
91: 
92: 

93 : 

94 : 
95: 

96 : 

97 : 
33 : 
9 9 : 

100 : 

101 : 
i n . 

X V x • 

^ r> — . 

x 'j _ : 

104 : 

105 : 

106 : 
107: 
103 : 
1 0 3 : 


112 ; 

i 7 •*: , 

. 

114 
Ilo ; 
116 
117 

■» *» O 

j. -i. o . 
11^ 
120 


123: 

124: 


399 


C 

c= 

c 


5 


IF( J.EQ. 999 ) STOP 

GO TO 379 

STOP 

END 
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SUBROUTINE GRAPHICS (J) 

DATA Y2M/0./ 

DO 1 1=1,201 
XI ( I) =T( I ) *1 .E+9 
Y2( I) =F3P( J.I) /1.E+? 

Y2M=AMAX1 C y5m,Y2 (I)) 

Yi(I)=FS( J.I) /l.E+9 

WRITE* 3.301) Y2M ~ 

FORMAT ( ' Max of graph=‘' ,G10 . 3 " Set it again for a ni 
+ Make it a little larger chan the original siax. T 

&eAd(3,) Y2M 

CALL MGRPHUXl'Y2*,loi?*Xl(l) . Xl(201),0., Y2M, &., 4.) 
CALL DATPL2 ( XI, Yl, 201,0^1,0,0. ,0,4,0) 

CALL TITLE C 1HL , 5Hl .E+9 , 5 , . 14 , 1 ) 

CALL TITLE ( 1HB,11HTIME ( nsec ) , 11 , . 14 , 1 ) 

CALL WHERE 

Z\jmU Jm l J XTi'i 

END ___ = __ = 

"sUBR5uriNE"MULTiGRA(NG)~~” , „ , , % 

COMMON / GRABLOCK / X ( 2 0 1 ) ,Y( 10,201) ,Z( 10,201) 

DIMENSION YG( 201 ) , YM(5) 

.5 scaling-— 

BO 3 J = 1 , NG 
DO 3 1=1,201 

YM< J)=AMAX1(YM< J) ,Y( J,I) ) 

DO 5 1=1,201 
DO 4 J = 1 , NG 
Y ( J , I ) = Y ( J , I ) / YM ( J ) 
x; I)=X( I)*1.E9 


CALL PINIT (1.5, 5., 1000) 


-PLOT 


YD ( T ) - Y 1 "* I ) 

CALL MGRP&L(X, YG.2Q1 ,X( 1 ) ,X(201) , 0., 1., 6, 
IF(NG.EQ.l) GO TO 999 


4. ) 


i-O • 

7-1 — . 

4. iL / • 


DO 6 
DO 3 

J = 2 , NG 
1=1,201 

123 : 

3 

YG( I ) 

=Y( J.I) 

129 : 

130 : 
1.5 1 : 
132: 

6 

CALL 

CALL 

CALL 

CALL 

DATPL2 C 
TITLE ( 1 
TITLE ( 1 
WHERE 

1x3 : 
134: 

939 

END 

IN 
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ce plot . " 
nax=?" ) 
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